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ABSTRACT 

Using data from the COSMOS survey, we perform the first joint analysis of galaxy-galaxy weak 
lensing, galaxy spatial clustering, and galaxy number densities. Carefully accounting for sample 
variance and for scatter between stellar and halo mass, we model all three observables simultaneously 
using a novel and self-consistent theoretical framework. Our results provide strong constraints on 
the shape and redshift evolution of the stcllar-to-halo mass relation (SHMR) from z = 0.2 to z = 1. 
At low stellar mass, we find that halo mass scales as Mh cx M° 46 and that this scaling does not 
evolve significantly with redshift from z = 0. 2 to z — 1. The slope of the SHMR rises sharply at 
M* > 5 x 10 10 M Q and as a consequence, the stellar mass of a central galaxy becomes a poor tracer 
of its parent halo mass. We show that the dark-to-stcllar ratio, M/j/M*, varies from low to high 
masses, reaching a minimum of M ft /M* ~ 27 at A/* = 4.5 x 10 10 M and M h = 1.2 x 10 12 M Q . 
This minimum is important for models of galaxy formation because it marks the mass at which the 
accumulated stellar growth of the central galaxy has been the most efficient. We describe the SHMR at 
this minimum in terms of the "pivot stellar mass", M* IV , the "pivot halo mass", M? lv , and the "pivot 
ratio", (Mh/M*) plv . Thanks to a homogeneous analysis of a single data set spanning a large redshift 
range, we report the first detection of mass downsizing trends for both Mf IV and M* IV . The pivot 
stellar mass decreases from Aff v = 5.75±0.13x 10 10 M at z = 0.88 to A/* piv = 3.55±0.17x 10 10 M Q at 
z — 0.37. Intriguingly, however, the corresponding evolution of A/ plv leaves the pivot ratio constant 
with redshift at (M^/M*) plv ~ 27. We use simple arguments to show how this result raises the 
possibility that star formation quenching may ultimately depend on M^/M* and not simply Mh, 
as is commonly assumed. We show that simple models with such a dependence naturally lead to 
downsizing in the sites of star formation. Finally, we discuss the implications of our results in the 
context of popular quenching models, including disk instabilities and AGN feedback. 
Subject headings: cosmology: observations - gravitational lensing - dark matter - large-scale structure 
of Universe, galaxies: evolution - stellar content 
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1. INTRODUCTION 

A fundamental goal in observational cosmology is to 
understand the link between the luminous properties 
of galaxies and the dark matter halos in which they 
reside. From an astrophysical perspective, measure- 
ments of the relationship between dark matter halo 
mass (Mh) and galaxy observables such as luminos- 
ity or stellar mass (M*) are critical for understanding 
how galaxy properties and their evolution with time 
are shaped by the halos that host them. Growing ev- 
idence suggests that halos accumulate stellar mass with 
an efficiency r](Mh,z) = {M„/Mh) x (flM/^b) that de- 
pends strongly on halo mass, peaking at Mh ~ 10 12 Ajf Q 
and dec lining towards lower and h i gher masses at z ~ 
(e.g.. iMandelbaum et all l2006bt IConrov k Wechsler 



2009; 



2010: 



Moster et a 
More et all 



.1120101: iBehroozT et al ,1 12010b IGuo et~a! 



20101 ). The stellar mass content is de- 



termined not only by the past merging of smaller sub- 
components but processes that, integrated over time, reg- 
ulate the conversion of gas into stars, including the rate 
at which fresh material is supplied to the halo, feedback 
mechanisms from supernovae (SN), galactic winds, and 
active galactic nuclei (AGN), and environmental effects 
such as ram pressure stripping, just to name a few. The 
global relationship between halo mass and average stel- 
lar content — the stellar-to-halo mass relation (SHMR) — 
probes the integrated outcome of these processes and 
as such, provides clues to their physical nature and 

Bower et all 



const r ains both semi-anal ytic models 
[200l iCroton et all l200l iSomerville eta! 
hydrodynamical simulations (iKeres et al 



2008TT and 
20051 I20QC 



Crain et al. 2009; Brooks et al.ll2009b iGabor et all 1201 



Agertz et al.N2011[) that aim to disentangle the relative 



contributions of such mechanisms 

From a cosmological perspective, the SHMR is vi- 
tal for determining how galaxies trace dark matter. 
A complete picture of the manner in which galaxies 
populate dark matter halos enhances the reconstruc- 
tion of the d ark matter power sp e ctrum from red- 
shift surveys (|Sanchez k Cold [2001 lYoo et all [2009]) 
and leads t o improved co n straints on cosmo l ogical 
parameters (lYoo et all 12006b iZheng k Weinberg! 12007b 
iCacciato et alJl2009f ). 

There are currently only two observational tech- 
niques capable of directly probing the dark matter 
halos of galaxies out to large radii (above 5 kpc) : 
galaxy-galaxy lensing (e.g.. iBrainerd et all IT996b 



12001b iHoekstra et al 



McKav et al 

2004b IMandelbaum et all 
2006ab Uohnston et alT ~ 



2006alht 



2004: Sheldon et al 



Hevmans et al 



Leauthaud et aL 20101) 



12007b ISheldon et all [2009b: 

and the kinematics of satel 



lite galaxies dMcKav et all 12002b iPrada et al 
Brainerd k Specianl 12003b van den Bosch et al 
W07J 



Conrov et all 1200 



Conrov 
2008b II 



More et all 



fBecker et al 
I2009L I2010D ~ 



12007b iNorberg et all 
The galaxy-galaxy 
g-g lensing") uses weak 



lensing technique (hereafter 
gravitational lensing to probe the gravitational potential 
around foreground ("lens") galaxies. The kinematic 

France-Hawaii Telescope with MegaPrime/MegaCam operated 
as a joint project by the CFHT Corporation, CEA/DAPNIA, the 
National Research Council of Canada, the Canadian Astronomy 
Data Centre, the Centre National de la Recherche Scientiflque 
de France, TERAPIX and the University of Hawaii. 



method uses satellite galaxies as test particles that trace 
the local velocity field (and hence the local gravitational 
potential). 

Another popular, albeit more indirect method, to in- 
fer the galaxy-dark matter connection is to measure 
the statistics of galaxy clustering. The results arc 
commonly interpreted using the halo occupation dis- 
tribution (HOD) model which describes the probabil- 
ity distribution P(N\Mh) that a halo of mass Mh is 
host to N galaxies abo ve some threshold in luminos- 
ity or stellar-mass (e.g.. iSeliald 120001: iPeacock k Smith! 
20MlScoccimarro et alll200lHBerlind k Weinberdl2002b 
Bullock et al.ll2002tlZehavi et aLll2002U2005HZheng et al l 
2005b 12007b iTinker et all 12007b IWake et alJ l201lF 
Zehavi et all 12010b iWhite et all l201lT Variations on 



the HOD approach include the conditional luminos- 
ity function §(L\Mh)dL which specifics the average 
number of galaxies of luminosity L zfc dL/2 that re- 



van den Bosch et al.ll2003b 


Vale k Ostrikeif2004blCooravl 


I'm Hi: Ivan den Bosch et al. 


120071 IVale & Ostrikerl 120081) 



and the conditional stellar mass function <$>(M*\M} l )dM* 
which yields the average number of galaxies with stellar 
masses in the rang e M* ± dM*/2 as a function of host 
halo mass M h (e.g.. [Yam ; et al.ll2009t iMoster et~aT] |2010l: 
iBehroozi et al.ll2010D . 

Finally, constraints on the SHMR have also been 
derived from the so-called "abundance matching" 
technique, which assumes that there is a monotonic cor- 
respondence between halo mass (or circul ar velocity) and 
galax y stellar mass (or l umino si ty) (e.g.. iKravtsov et all 



2001 IVale fc Ostrikerl 120041: iTasitsiomi et all 



Vale fe Ostrikerl 12006b IConrov k Wechsler 
Drorv et all 12009b IMoster et all 12010b IBehroozi et all 



2010b IGuo et alll2010D ~ 



While individual applications of the techniques de- 
scribed above have p rovided important insight, in 
iLeauthaud et al.l ([201 ID (hereafter "Paper I") we take 
a further step by combining separate probes into a 
self-consistent theoretical framework. Specifically, our 
method combines measurements of g-g lensing, galaxy 
clustering, and the galaxy stellar mass function. Begin- 
ning with the standard HOD formalism, we make several 
modifications that enable us to a) extract the parameters 
that determine the SHMR and to b) simultaneously fit 
data from multiple probes while allowing for independent 
binning schemes for each probe. The goal of the current 
paper is to apply the methodology developed in Paper I 
to observations from the COSMOS survey from z = 0.2 
to z = 1.0. This enables confident measurements of the 
shape of the SHMR and its evolution with time, with 
important implications for models of galaxy formation. 

While stellar mass estimates are a key galaxy observ- 
able in this work, it is important to highlight the un- 
certainties and sometimes unknown systematic biases 
that affect them and, if not treated carefully, can mud- 
dle attempts to c ompare results from disparate surveys 
(see discussion in IBehroozi et al1l2010l ). One of the ad- 
vantages of the COSMOS data set is that evolution- 
ary trends can be studied within the sample using self- 
consistent stellar mass estimates, making our conclusions 
more robust. 

The layout of this paper is as follows. The data are 
described in $2] followed by the presentation of the g-g 
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lcnsing, clustering, and SMF measurements in In the 
interest of brevity, we only give a short and necessarily 
incomplete review of the theoretical background in S{4j 
We strongly encourage the reader to refer to Paper I for 
a complete description of the theoretical foundations for 
this work. Our main results are presented in iJU Finally, 
we discuss the results and draw up our conclusions in 8E1 
and 83 

We assume a WMAP5 ACDM cosmology with fi m = 
0.258, fl A = 0.742, fl h h 2 = 0.02 273, n s = 0.963 <t 8 = 
0.796, H = 72 km s" 1 Mpc" 1 (Iffinshaw et all 120091) . 
All distances are expressed in physical Mpc units. The 
letter Mh denotes halo mass in general whereas M2oob is 
explicitly defined as M 2 oob = M{< r 2 oob) = 200p|7rr| oofc 
where r2oob is the radius at which the mean interior den- 
sity is equal to 200 times the mean matter density (/?). 
Stellar mass is noted M* and has been derived using a 
Chabrier Initial Mass Function (IMF). Stellar mass scales 
as 1/Hq. Halo mass scales as 1/H . All magnitudes are 
given on the AB system. 

2. DESCRIPTION OF DATA 

The COSMOS survey (IScoville et all 12001 brings to- 
gether a broad array of panchromatic observations with 
imaging data from X-ray to radio wavelengths and a 
large spect roscopic follow-up program (zCOSMOS) with 
the VLT ({Lilly et alJl200l . In particular, the COS- 
MOS program has imaged the largest contiguous area 
(1.64 degrees 2 ) to date with the Hubble Space Telescope 
(HST) using the Advanced Camera for Surveys (A CS) 
Wide Field Channel (WFC) (jKoekemoer et alil2007D . 



2.1. The ACS Tensing catalog 

The general methodology for the construction of the 
COSMOS ACS weak lensing catalog an d our shape mea- 
surement p rocedure are p resented in ILeauthaud et all 
(|2007D and iRhodes et aTT (|2007l) . In this section, we 
present several updates to the pipeline that we have im- 
lemented since those publications. In ILeauthaud et all 
2001 we used a parametric correction for the effects 
of Charge Transfer Inefficiency (CTI) on galaxy shape 
measurements. Instead, in this paper, we use a physi- 
cally motivated CTI correction scheme that operates on 
the raw data and returns electrons to pixels from which 
they were unintentionally dragged out d uring readout. 
This c orrection scheme has been shown by iMassev et all 
(|2010T ) to reduce the CTI trails by a factor of -10 every- 
where in the CCD and at all flux levels. 

Following CTI correction in the raw images, image reg- 
istration, geometric distortion, sky subtraction, cosmic 
ray rejection and the final combination of the dithered 
images are performed b y the MultiDriz zle algorithm 
Koekemoer et al1l2002f ). As described in IRhodes et all 



20071) . a finer pixel scale of 0.03"/pix was used for 
the final co-added images. The lensing source cata- 
log is constructed from 575 ACS/WFC tiles. Defects 
and diffraction spikes are carefully removed, leaving a 
total of 1.2 x 10 6 objects to a limiting magnitude of 
If814W = 26.5. 

The next step is to measure the shapes of galaxies 
and to correct them for the distortion induced by the 
time varying ACS PSF (Point Spread Function) (see 
IRhodes etldl 120071) . We continue to use a PSF model 



based on physical par ameters rather than arbi trary prin- 
cipal components. In ILeauthaud et a l. (2007), we mod- 
elled the PSF as a multivariate polynomial in x, y and 
focus. We now fit the PSF as a function of x, y, fo- 
cus and velocity aberration of the poi nting (recorded in 
CALA CS headers as "VAFACTOR"). ISchrabback et all 
(|2010D found that VAFACTOR partially correlates with 
the higher-order PSF variations, motivating our use of 
this quantity. In the g-g lensing analysis presented here, 
the weak lensing shear is azimuthally averaged. Thus, 
any effects of PSF anisotropy cancel to leading order. 
Therefore, our science analysis is insensitive to subtle 
differences in the PSF modeling. We confirmed this by 
repeating our analysis with the independently o btained 
weak lensing catalog by ISchrabback et all (|2010f ) , which 
yields fully consistent results. 

Finally, simulated images are used to derive the shear 
susceptibility factors that are necessary in order to trans- 
form shape measuremen ts into unbiased shear estimators 
({Leauthaud et al.ll2007l) . Representing a number density 
of 66 galaxies per arc-minute 2 , the final COSMOS weak 
lensing catalog contains 3.9 x 10 5 galaxies with accurate 
shape measurements. 

2.2. Photometric and Spectroscopic Redshifts 

We use two updated versions (vl.8 dated from the 13 th 
of July 2010 and vl.7 dated from the 1 st of August 2009) 
of the p hotometric r e dshift s (hereafter photo-z's) pre- 
sented in lllbert et al.l (|2009l ) which have been computed 
with over 30 bands of multi-wavelength data. In par- 
ticular, deep K s . J. and u* band data allow for a good 
photo- z estimate at z > 1 via the 4000A break which is 
increasingly shifted into the near infra-red (IR). Further 
details r egarding the data and the photometry can be 
found in lCapak et""aTl (|2007l ). 

The photo-z catalog vl.8 has improved redshifts at z > 
1 compared to vl.7. At z < 1, the difference between 
the two catalogs is minor. We use catalog vl.8 for g-g 
lensing measurements and vl.7 for the SMF and galaxy 
clustering. Interchanging the two catalogs does not affect 
our results. 

Photo-z 's were estimated using a \ 2 template fit- 
ting method (Le Phare) and calibrate d with large spec- 
troscopic sample s from VLT-VIMOS ({Lilly et al.l 120071: 
lllbert et al.ll2009D and Keck-DEIMOS. The dispersion in 
the photo-z 's as measured by comparing to the spectro- 
scopic redshifts is o"A z /(i+ Za ) = 0.007 at i\ B < 22.5 
where Az = z spcc — z p \ wt . The deep IR and IRAC (In- 
frared Array Camera on the Spitzer Space Telescope) 
data enables the photo-z's to be calculated even at fainter 
magnitudes with a reasonable accuracy of CA2/(i+z Bpoc ) = 
0.06 at ito ~ 24. 

Figure |T] compares the spectroscopic and photometric 
redshifts of 8812 galaxies that belong to both the lens- 
ing source catalog and the zCOSMOS "bright" or "faint" 
programs. This Figure also illustrates the sensitivity of 
g-g lensing signals to photometric redshift errors. As can 
be seen from Figure [IJ g-g lensing signals are increas- 
ingly insensitive to photometric redshift errors for source 
galaxies at higher redshifts. 

2.3. Stellar Mass Estimates 
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Fig. 1. — The effects of photometric rcdshifts errors on galaxy- 
galaxy lcnsing signals. This figure illustrates the quality of the 
photometric rcdshifts for source galaxies in this paper by compar- 
ing them to a combined sample of 8812 spectroscopic rcdshifts from 
the zCOSMOS "bright" and "faint" programs for which we have 
applied the same selection as for source galaxies. We emphasize 
that since the spectroscopic sample does not go as faint as our 
source catalog, this figure is necessarily an underestimate of the 
true rcdshift errors for the source catalog - a complete sample of 
faint galaxies with spectroscopic rcdshift would be necessary in or- 
der to fully test the photometric rcdshifts. There are three ways in 
which photometric rcdshift errors can impact g-g lensing signals. 
Firstly, any type of photometric error such that z p hot < z\ ena will 
have no effect on the signal because such objects are not included 
in the background selection (bottom hashed region). Secondly, 
photometric errors such that z p hot > zi ena and z spe c < ziens wm 
lead to a signal dilution (left hashed region). Finally, photometric 
errors such that z phot > z lcns and z apcc > z lcns but z phot ^ z spcc 
will lead to a bias in AS because S cr ; t will be mis-estimated when 
transforming 7 into AS. The dotted and dashed lines indicate 
where photo-z errors lead to a 10% and a 20% error on AE for 
a lens located at z=0.2. As can be seen, g-g lcnsing signals are 
increasingly insensitive to photometric rcdshift errors for source 
galaxies at higher redshifts. 

Stellar ma sses are estim a ted u sing the Bayesian code 
described in iBundv et al.1 (|2006t ) assuming a Chabrier 
IMF. Briefly, an observed galaxy's spectral energy dis- 
tribution (SED) and redshift is referenced to a grid o f 
models constructed using the lBruzual fc Chariot! pOOl 
synthesis code. The grid includes models that vary in 
age, star formation history, dust content, and metallic - 
ity The assumed dust model is iCharlot fc Fall (|2000( ). 
At each grid point, the probability that the observed SED 
fits the model is calculated, and the corresponding stel- 
lar mass to K-band luminosity ratio and stellar mass is 
stored. By marginalizing over all parameters in the grid, 
the stellar mass probability distribution is obtained. The 
median of this distribution is taken as the stellar mass 
estimate, and the width encodes the uncertainty due to 
degeneracies and uncertainties in the model parameter 
space (also described as "model error" in § 14. 2[) . The 
final uncertainty on the stellar mass also includes the 



K-band photometry uncertainty as well as the expected 
error on the luminosity distance that results from the 
photo-z uncertainty. The typical final uncertainty is 0.1- 
0.2 dex (also see § 14.21 and Figure [4|). Systematic uncer- 
tainties also come from the choice of stellar population 
templates used to fit the data as well as the assumption 
of a universal IMF. To first order these will cause global 
offsets in the mass estimates that will not affect compar- 
isons within our sample but may impact comparisons to 
work by other authors. Since the primary goal of this 
paper is to study the redshift evolution of the SHMR 
derived from COSMOS data alone, we do not include 
th ese systematic uncer tainti es in our analy s is and refer 
to iConrov et all (|2009f ) and iBehroozi et~afl (|2010l) for a 
broad discussion of systematic errors in stellar mass es- 
timates. 

Following the approach in IBundv et al.l (j2010l ) , we ob- 
tained PSF-matched 3'.'0 diameter aperture photome- 
try from the ground-based COSMOS cat alogs (filters 
u*,B,j , y,j,g + ,r + , i + , z + , K s ) described in ICapak et al.1 
(|200l . Inbert^TaLl (I2009I) . and lMcCracken et al.l (12010^ 7 
after ap plying the ph o tomet ric zero-point offsets tabu- 
lated in ICapak et al.l (|2007f ). The depth in all bands 
reaches at least 25th magnit ude (AB) with the JT. 9 -band 
limited to K s < 24. Unlike iDrorv et al.l (|2009f ). we re- 
quire if s -band detections for all galaxies i n the sample. 
We ha ve found that the mass estimates in | Bundv et al.l 
(|201Q[) agree with those of lDrorv et~aTl (|2009f ) within the 
expected uncertainties (i.e., < 0.2 dex). The mass esti- 
m ates used in this wo rk arc slightly different from those 
in lBundv et alj (|2010f ) in that they are based on updated 
rcdshift information (vl.7 of the photo-z catalog and the 
latest available spectroscopic redshifts as compiled by the 
COSMOS team) and use a slightly different cosmology 
(H = 72 km s" 1 Mpc" 1 instead of H Q = 70 km s" 1 
Mpc" 1 ). 

The bins in our analysis are defined by using two esti- 
mates of the mass completeness of the sample, as deter- 
mined by the magnitude limits K s < 24 and isi4W < 25. 
The first estimate is more conservative and comes from 
estimating the observed magnitude of a maximal M*/L 
stellar population model with solar metallicity, no dust, 
and a r = 0.5 Gyr burst of star formation that occurred 
at Zform = 5. As a function of rcdshift, the stellar mass of 
such a population at the point where its observed K s - and 
-^8i4W-t>and flux falls below the magnitude limits defines 
the mass completeness and roughly matches the 80% 
completen ess limits determine d when deeper samples are 
available (|Bundv et al.l l2006h . This redshift-dependent 
limit is plotted as the dashed line in Figure [5J In prac- 
tice, low mass galaxies exhibit more star formation and 
therefore have lower M*/L ratios than the passive tem- 
plate described above so we also define a second, more 
liberal mass limit (i.e., lower) that corresponds to a star- 
forming population. This is plotted as the dotted line. 
The majority of the sample bins lie above the more con- 
servative limit, but the lowest mass bins are allowed to 
reach the star-forming mass limit under the assumption 
that passive stellar populations (potentially missed) at 
such low masses are extremely rare. 

3. MEASUREMENTS 
3.1. Sample Selection 
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We use two COSMOS galaxy catalogs: the Subaru 
(photo-z) catalog and the ACS (lensing) catalog. The 
ACS catalog corresponds to the COSMOS area that has 
been imaged with HST and covers a slightly smaller area 
than the Subaru catalog. The ACS coverage of COSMOS 
is 1.64 degrees 2 and the Subaru coverage of COSMOS is 
2.3 degrees 2 . The ACS catalog is used for our g-g lensing 
measurements and for the SMFs. The Subaru catalog is 
used to calculate the galaxy clustering. 

Galaxies arc selected with K s < 24 to ensure that stel- 
lar masses can be computed for all galaxies in the sam- 
ple. Note that for the g-g lensing analysis, these cuts 
only apply to the foreground lens sample. In addition 
to these cuts, we also reject galaxies that are in masked 
areas where the photometry is deemed unreliable. For 
this, we use the union of four masks in the B, V, ip, and 
zp bands ("COSMOS.B.mask", "COSMOS.V.mask" , 
"COSMOS.ip.mask", "COSMOS.zp.mask"). 

For the ACS sample, star-galaxy separation is per- 
formed using th e morp hological classifier described in 
iLeauthaud et al.l (|2007f ). For the Subaru sample, we re- 
ject stars by imposing the condition that Xg a i — Xstar > 
0.5 where x 2 al represents the chi-square value of the best- 
fitting galaxy SED template and x 2 tar represents the chi- 
square value of the best-fitting stellar SED template. In 
conjunction with the K s cut, this stellar selection agrees 
with our ACS based morphological classification at the 
97% level. According to our ACS morphological classi- 
fier: 2% of the objects in the final Subaru catalog are 
classified as stars and 0.1% of galaxies are mis-classified 
as stars. 

We compute the g-g lensing, galaxy clustering, and 
stellar mass in three redshift intervals: z\ = [0.22,0.48], 
z 2 = [0.48,0.74], and z 3 = [0.74,1]. For each redshift 
bin, we define a lower limit on the stellar mass to ensure 
that all samples are complete in terms of stellar mass. 
These cuts correspond to M* > 10 8 - 7 M Q for z 1} M* > 
10 9 - 3 Mq for z 2 , and M* > 10 9 8 M Q for z 3 (see Table [J) 
and are defined by the stellar mass completeness limit at 
the far edge of the redshift bin (see Figured]). All galaxy 
samples used in this paper are complete. 

3.2. Stellar Mass Functions 

The SMFs are calc ulated using the ACS catalog. As 
in lBundv et al.l (12010 1 ) . we compute mass functions using 
the V m ax technique |Schmidtlll968D . We weight galaxies 
by the maximum volume in which they would be de- 
tected within the K s -band and Isi4w-band limits in a 
given redshift interval. For each host galaxy i in the 
redshift interval j, the value of V 7 l nax is given by the min- 
imum redshift at which the galaxy would drop out of the 
sample, 



where the redshift interval, j, is defined by \z 3 min , z 3 maa \ 



and 



V 



dz 



(1) 



2io 



where dQ is the solid angle subtended by the survey area, 
and dV/dz is the comoving volume element. The redshift 
limits are given as, 



Zhigh — ninfioi! Z K lim ' Z I U 



(2) 
(3) 



and 



refer to the redshift at which the 



galaxy would still be detected below the K s - and Isi4W- 
band limits. We use the best-fit SED template as deter- 
mined by the stellar mass estimator to calculate these 
values, thereby accounting for the fc-corrections neces- 
sary to compute V max values (no evolutionary correction 
is applied). 

3.3. Galaxy Autocorrelation Function 

The galaxy clustering samples are defined by a series of 
stellar mass thresholds rather than bins. For galaxy clus- 
tering, threshold samples are the most straightforward to 
model in the HOD context. At small scales, the signal-to- 
noise (S/N) in galaxy clustering is due mostly to Poisson 
noise and at large scales it is subject to sample ( "cosmic" ) 
variance. The optimal binning scheme for clustering can 
be different than that for g-g lensing. For example, mas- 
sive galaxies produce the strongest shear signals, thus a 
relatively small sample is required to produce a robust 
measurements relative to a simple pair-counting statis- 
tic like the two-point correlation function whereas the 
inverse would be true for low-mass galaxies. The thresh- 
old samples we employ for clustering measurements are 
listed in Table [3J Except at the high mass end where 
galaxies become scarce, the binning scheme employed for 
the clustering is a constant 0.5 dex in log 10 (M*). 

Since we do not need galaxy shape information to mea- 
sure the angular correlation function w(9), we are not re- 
stricted to the ACS coverage of the COSMOS field. Wc 
therefore use the COSMOS Subaru catalog to calculate 
w(9). For each threshold sample we measure w(9) using 
the well-known Landy-Szalay estimator of 



DD - 2DR + RR 
RR ' 



(4) 



where DD are the number of data-data pairs in a given 
bin of angle 6, RR are the number of random-random 
pairs, and PR are th e number of data-random pairs 
(jLandv fc Szalavlll993| ) . Data and random pairs are nor- 
malized by the total number of galaxies and randoms, 
respectively. In all measurements we use 10 5 randoms. 
The distribution of the randoms are taken from a com- 
bination of four COSMOS masks (B, V, ip, zp), thus 
mimicking the angular completeness and geometry of the 
survey. 

Because the volumes probed in each sample are small, 
our clustering measurements a re subject to the effec t 
of the integral constraint (IC) (jGroth fc Peebleslll977ft . 
Due to spatial fluctuations in the number density of 
galaxies, the mean correlation function measured from an 
ensemble of samples will be smaller than the correlation 
function measured from a single contiguous sample of the 
same volume as the sum of the ensemble sample. This 
attenuation of w(8) becomes relevant on angular scales 
significant with respect to the sample size. We estimate 
the IC correction to our w(9) measurements through the 
use of mock galaxy distributions described in Paper I and 
§21 In practice, we do not modify the measurements for 
the IC but rather adjust the theoretical models to ac- 
count for the finite sample size. 

Finally, in order to compute w(9) from our model, we 
need to know the normalized redshift distribution of the 
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TABLE 1 

Characteristics of three redshift bins 



Zmin Zmax z median ACS Volume Subaru Volume Min M* N galaxies ACS a N galaxies Subaru a 
10 6 Mpc 3 10 6 Mpc 3 



0.22 0.48 0.37 0.88 1.24 10 8 7 M Q 14956 20426 

0.48 0.74 0.66 2.03 2.84 10 9 3 M© 15103 20068 

0.74 1.0 0.88 3.13 4.39 10 9 8 M© 14387 18853 



a Numbers are quoted after stars have been removed, the selection K s < 24 as well as the lower limit stellar mass 
limit has been applied, and objects in masked areas have been removed from the catalog. 
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0.0 0.2 0.4 0.6 0.8 1.0 

Photometric Redshift 

Fig. 2. — Distribution of COSMOS galaxies as a function of redshift and stella r ma ss. The completeness limit for passive galaxies is 
shown by the red dashed line and for active galaxies by the blue dotted line (sec § 12.31 for the derivation of the completeness limits) . The 
green hashed region approximately shows the cross-over mass above which the number of red galaxies becomes larger than the number of 
blue galaxies. The vertical solid dark blue lines delineate the three redshift bins Zi, 221 and 23 while the horizontal solid dark blue lines 
show the binning scheme for the g-g lensing measurements. 



TABLE 2 

Thresholds for the angular correlation function 







w binl 


w bin2 


w bin3 


w bin4 


w bin5 


w bin6 


minimum log 10 (M„), z\ = 


0.22,0.48] 


11.1 


10.8 


10.3 


9.8 


9.3 


8.8 


minimum log 10 (M»), 22 = 


0.48,0.74] 


11.1 


10.8 


10.3 


9.8 


9.3 


none 


minimum log 10 (A4"*), 23 = 


0.74, 1.0] 


11.1 


10.8 


10.3 


9.8 


none 


none 



galaxy sample, N(z) (see Equation 32 in Paper I). For 
this we use the probability distribution functions of the 
photometric redshifts to estimate the true N(z) for our 
photo-z slices. However, the COSMOS photo-z errors 
are accurate enough such that in tests we find a minimal 
effect when using a flat top-hat z-bin for N(z). 

3.4. Galaxy-galaxy lensing: from galaxy shapes to AE 

In the weak gravitational lensing limit, the observed 
shape e bs of a source galaxy is directly related to the 
lensing induced shear 7 according to 



£obs = £int + 7, (5) 

where £i n t is the source galaxy's intrinsic shape that 
would be observed in the absence of gravitational lens- 
ing. In our notation, £j nt , £ bs, and 7 are spin-2 tensors. 
The above relationship indicates that galaxies would be 
ideal tracers of the distortions caused by gravitational 
lensing if the intrinsic shape £i n t of each source galaxy 
was known a priori. However, lensing measurements 
exhibit an intrinsic limitation, encoded in the width 
of the ellipticity distribution of the galaxy population, 
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noted here as ai n t, and often referred to as the "intrinsic 
shape noise " . Because the intrinsic shape noise (of order 
Omt ~ 0.27. lLeauthaud et al.ll2007t) is significantly larger 
than 7, shears must be estimated by averaging over a 
large number of source galaxies. 

Throughout this paper, the gravitational shear is noted 
as 7 whereas 7 represents our estimator of 7. The un- 
certainty in the shear estimator is a combination of un- 
avoidable intrinsic shape noise, cr^t = yf (sf nt ), and of 
shape measurement error, <7 mea s: 



(6) 



We will refer to as "shape noise" whereas a ln t will 
be called the "intrinsic shape noise" . The former includes 
shape measurement error and will vary according to each 
specific data-set and shape measurement method. Aver- 
aged over the whole COSMOS field, the weak lensing 
distortions represent a negligible perturbation to Equa- 
tion [HI 

The derivation of our shear estimator is presented in 
ILeauthaud et al.l (120071). We employ the RRG method 
(see iRhodes et aL 120001 for further details) for galaxy 
shape measurements. Briefly, we form 7 from the PSF 
corrected ellipticity according to 



C x 



^obs 

~G 



(7) 



where the shear susceptibility factoiQ, G, is measured 
from moments of the global distribution of e b s and 
ot her, higher order sh ape parameters (see Equation 28 
in IRhodes et al.l 120001 ) . Using a set of simulated im- 
ages similar to those of Shear TEstin g Pro gram (STEP; 
iHevmans et al.l l2006ri iMassev et all 120071) but tailored 
exclusively to this data-set, we find that, in order to cor- 
rectly measure the input shear on COSMOS-like images, 
the RRG method requires an overall calibration factor of 

G = (0.86+gS)- 1 - . 
The shear signal induced by a given foreground mass 

distribution on a background source galaxy will depend 
on the transverse proper distance between the lens and 
the source and on the redshift configuration of the lens- 
source system. A lens with a projected surface mass 
density, £(r), will create a shear that is proportional to 
the surface mass density contrast, AE(r): 



AE(r) = £(< r) - S(r) = S crit x 7t (r). 



(8) 



Here, £(< r) is the mean surface density within proper 
radius r, S(r) is th e azimuthally averaged surface den- 
sity a t radius r (e.g, lMiralda-Escudelll991l ; iWilson et al.l 
120011 ). and 74 is the tangentially projected shear. The ge- 
ometry of the lens-source system intervenes through the 
critical surface mass density S C rit, which depends on 
the angular diameter distances to the lens (£>ol), to the 
source (Dos), and between the lens and source (£>ls) : 



^os 



(9) 



4ttGn D ol D ls 

1 Not to be confused with Newton's constant which we have 
noted Gn. 

2 Note that some authors define the comoving critical surface 
mass density which has an extra factor of (1 + z)~ 2 with respect 
to ours due to the use of comoving instead of physical distances. 



where Gn represents Newton's constant. Hence, if red- 
shift information is available for every lens-source pair, 
each estimate of jt can be directly converted to an esti- 
mate of AE(r) which is a more desirable quantity than 
7t because it depends only on the mass distribution of 
the lens. 

To measure AE(r) with high signal-to-noise, the lens- 
ing signal must be stacked over many foreground lenses 
and background sources. For every zth lens and jth 
source separated by a proper distance ry, an estimator 
of the mean excess projected surface mass density at r.y 
is computed according to: 



1t,ij ^ ^crit.zj? 



(10) 



where 7^ is the tangential shear of the source rela- 
tive to the lens. The COSMOS photometric redshifts 
described in W2.'2\ are used to estimate £ c rit,ij f° r every 
lens-source pair. In order to optimize the signal-to-noise, 
an inverse variance weighting scheme is employed when 
AEjj is summed over many lens-source pairs. Each lens- 
source pair is attributed a weight that is equal to the 
estimated variance of the measurement: 



1 



(Ecrit,i 



iY 



(ii) 



In this manner, faint small galaxies which have large 
measurement errors are down-weighted with respect to 
sources that have well measured shapes. 

For the types of lenses studied in this paper, the signal- 
to-noise per lens is not high enough to measure AT, on 
an object by object basis so instead we stack the signal 
over many lenses. For a given sample of lenses, the total 
excess projected surface mass density is the weighted sum 
over all lens-source pairs: 



AS 



7=1 2^i=l W ij x lt,ij X Sorit.i 



3.5. Galaxy-Galaxy Lensing Measurements 



N S a 



(12) 



We only give a brief outline of the overall methodol- 
ogy used to compute the g-g lensing signals since this 
has al ready been presented in detail in ILeauthaud et al.l 
(|2010D . Foreground lens galaxies are divided into three 
redshift samples and then arc further binned by stellar 
mass (see Figure and Table For each lens sample, 
AE is computed according to Equation [T^] from 25 kpc 
(physical distance) to 1. 5 Mpc in logarithmicall y spaced 
radial bins of 1.8 dex. In lLeauthaud et ail (|2010l ) we used 
a theoretical estimate of the shape measurement error 
in order to derive the inverse variance for each source 
galaxy. Instead, in this paper, the dispersion of each 
shear component is measured directly from the data in 
bins of S/N and magnitude. The measured shear dis- 
persion is equal to the quadratic sum of the intrinsic 
shape noise and of the shape measurement error. Our 
empirical derivation of the shear dispersion varies from 
cr~ ~ 0.25 for bright galaxies with high S/N to ~ 0.4 
for faint galaxies with low S/N. Overall, we find that 
the theoretical and the empirical schemes yield very sim- 
ilar results with the latter method resulting in slightly 
larger error bars because the theoretical scheme tends to 
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somewhat underestimate the shape measurement error 
for faint galaxies. 

Photometric redshifts are used to derive £ C rit for each 
lens-source pair. The lower 68% confidence bound on 
each source redshift is used to select background galaxies. 
For each lens-source pair, we demand that Source - 2i e ns > 

source 

) so as to minimize foreground contamina- 
tion. The g-g lensing signal is most sensitive to redshift 
errors when z SOU rce is only slightly larger then zi ons (see 
Figure [lj . For this reason, in addition to the previous 
cut, we also implement a fixed cut so that z SO urcc — zicns > 
0.1. Furthermore, in order to minimize the effects of sig- 
nal dilution caused by catastrophic errors, we also reject 
all source galaxies with a secondary peak in the redshift 
probability distri bution function (i. e. the parameter ZP2 
is non zero in the lllbert et al.ll2009l catalog). This cut is 
aimed to reduce the number of catastrophic errors in the 
source catalog. After all cuts have been applied, the g-g 
lensing source catalog that we use represents 35 galaxies 
per arcmin 2 . 

Final ly, we re-compute all ou r g-g lensing signals us- 
ing the lSchrabback et all (|2010f ) COSMOS shear catalog 
which has been independently derived from ours. We find 
identical g-g lensing signals from both shear catalogs, in- 
dicating that any relative shear calibration differences 
between the two shear catalogs has no impact on these 
results. This test provides an independent validation of 
our g-g lensing results. 

4. THEORETICAL FRAMEWORK 

Paper I presents the general theoretical foundations 
that form the backbone of this paper. In this section, we 
only give a brief, and thus necessarily incomplete, review 
of the theoretical background and strongly encourage the 
reader to refer to Paper I for further details. We adopt 
the same model and notation as in Paper I. 

Paper I describes an HOD-based model that can be 
used to analytically predict the SMF, g-g lensing, and 
clustering signals. The key component of this model is 
the SHMR which is modelled as a log-normal probability 
distribution function with a log-normal scatteiQ noted 
ciogM, an d with a mean-log relation denoted as M* = 
/sHMa(Af fc ). 

For a given parameter set and cosmology, / SH mr and 
ciogM, can be used to determine the central and satellite 
occupations functions, (N cen ) and (iV sa t). These are used 
in turn to predict the SMF, g-g lensing, and clustering 
signals. 

4.1. The stellar-to-halo mass relation 

Following iBehroozi et~aT1 (|2010f ) (hereafter "BIO"), 
/shmr(-^/i) is mathematically defined via its inverse func- 
tion: 

logio(/^La(M„)) = Iog 10 (M h ) = 

(13) 

3 Scatter is quoted as the standard deviation of the logarithm 
base 10 of the stellar mass at fixed halo mass. 



where M\ is a characteristic halo mass, M* t o is a charac- 
teristic stellar mass, j3 is the faint end slope, and 6 and 7 
control the massive end slope. We refer to BIO for a more 
detailed justification of this functional form. Briefly, we 
expect that at least 4 parameters are required to model 
the SHMR: a normalization, break, a faint end slope, and 
a bright end slope. In addition, BIO have found that the 
SHMR displays sub-exponential behaviour at large M*. 
This is modelled by the 8 parameter which leads to a to- 
tal of 5 parameters. Figure [3] illustrates the influence of 
each parameter on the shape on the SHMR and further 
details on the role of each parameter can be found in § 
2.1 of Paper I. 

In contrast to BIO, we do not parameterize the red- 
shift evolution of this functional form. Rather, we bin the 
data into three redshift bins and check for redshift evolu- 
tion in the parameters a posteriori. We also assume that 
Equation [13] is only relevant for central galaxies. Follow- 
ing the HOD ansatz, satellite galaxies within groups and 
clusters occupy subhalos — bound, virialized halos that 
are contained within the radius of a larger halo. Abun- 
dance matching models like BIO assume that the halos 
and subhalos of the same mass (or circular velocity) con- 
tain galaxies of the same stellar mass, where the subhalo 
mass is taken as the mass at the last time the galaxy 
was a central galaxy. Here we put no such prior on the 
galaxy-subhalo connection. Rather, the halo occupation 
of satellite galaxies is constrained by the data. 

4.2. Scatter between stellar and halo mass 

The measured scatter in stellar mass at fixed halo mass 
has an intrinsic component (noted crj M ^), but also in- 
cludes a stellar mass measurement error due to redshift, 
photometry, and modeling uncertainties (noted c{£ Mji ). 
Ideally, we would measure both components but unfortu- 
nately we can only constrain the quadratic sum of these 
two sources of scatter. Nonetheless, given a model for 
°iogM, ' we cou ld m principle extract c[ ogMt from ai og M, • 

Previous work suggests t hat owm. i s inde pendent of 
halo mass. For exampl e. lYang et al.l (I2009T ) find that 
"logM, = 0-17 dex and iMore et alJ (|2009l )find a scat- 
ter in luminosity a t fixed halo mass of 0.16 ± 0.04 dex. 
Both iMoster et all (|2010l ) and B10 are able to fit the 
SDSS (Sloan Digital Sky Survey) galaxy SMF assum- 
ing <7io g M» = 0.15 dex and cria S M» = 0.175 dex respec- 
tively. However, these results are derived with spectro- 
scopic samples of galaxies. In contrast to these surveys, 
we expect a larger measurement error for the COSMOS 
stellar masses due to the use of photometric redshifts. In 
addition, since photo- z errors increase for fainter galax- 
ies, we might also expect that c™ gMi> (and thus eri og M„) 
will depend on M*. 

To test if the assumption that <7io g M, is constant has 
any impact on our results, we implement two models for 
ciogM, • I n the first case (called "sig_mod1"), erio g M„ is 
assumed to be constant (this is our base-line model). In 
the second case (called "sig_mod2"), wc explicitly model 
c^iogM to reflect stellar mass measurement errors. Note 
that the goal of this exercise is not to perform a careful 
and thorough error analysis, but simply to build a realis- 
tic enough model to asses whether or not a M* dependent 
error has any strong impact on our conclusions. 

For the SIG_mod2 model, we consider three contribu- 
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TABLE 3 
Binning scheme for the g-g lensing 











g-g binl 


g-g bin2 


g-g bin3 


g-g bin4 


g-g bin5 


g-g bin6 


g-g bin7 


*1 


= [0.22,0.48] 


min 

max 


Iog 10 (Af„) 
log 10 (Af„) 


11.12 
12.0 


10.89 
11.12 


10.64 
10.89 


10.3 
10.64 


9.82 
10.3 


9.2 
9.8 


8.7 
9.2 


Z2 


= [0.48.0.74J 


nun 

max 


log 10 (M.) 
log 10 (Af„) 


11.29 
12.0 


11.05 
11.29 


10.88 
11.05 


10.65 
10.88 


10.3 
10.65 


9.8 
10.3 


9.3 
9.8 


Z3 


= [0.74,1.0] 


mm 
max 


Iog 10 (Af*) 
log 10 (M*) 


11.35 
12.0 


11.16 
11.35 


10.97 
11.16 


10.74 
10.97 


10.39 
10.74 


9.8 
10.39 


none 
none 
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Fig. 3. — Illustration of the influence of Mi, M* } o, 0, S and 7 on the shape of the SHMR. 



tions to the stellar mass error budget. The first is called 
"model error" : this is measured by the 68% confidence 
interval of the mass probability distribution determined 
for each galaxy by the mass estimator. It represents the 
range of model templates (each with its own M/L ratio) 
that provide reasonable fits to the observed SED. This 
range is determined by the photometric uncertainty in 
the observed SED, degeneracies in the grid of models 
used to fit the data, and how well the grid of models 
represents the true parameter space of observed galaxy 
populations as well as their colorsfj The second term 
is the photo- 2: error, which derives from the uncertainty 
in the luminosity distance owing to the error on a given 
photometric rcdshift. The final component is the pho- 
tometric uncertainty from the observed K-band magni- 
tude, which translates into an uncertainty in luminosity 
and therefore stellar mass. The total measurement error, 
^logM, i s t ne sum m quadrature of these three sources of 
error. The results are shown in Figure 2] for the three 
rcdshift bins. 

As detailed in § [5j we find that our results are largely 
unchanged, regardless of which form we adopt for CTi og M, ■ 
This can be explained as follows. Since the data are 
binned by M* , the observables are in fact sensitive to the 
scatter in halo mass at fixed stellar mass, noted cn og M h ■ 
Given a model for the SHMR, <Ti og M h can be mathe- 
matically derived from cri og M, ■ Further details on the 

4 Model template problems, especially in the restframe near-IR, 
and agc-mctallicity degeneracies may help explain the slight rise in 
the "model error" contribution to the total stellar mass uncertainty 
for high mass galaxies at low redshift in Figure 4. 



mathematical connection between between cri og M, and 
cio g M h can be found in Paper I. We find that the slope 
of the SHMR increases steeply at A/* > 10 11 M Q so 
that cri og M h becomes quite large at the high mass end. 
For example, cri og M h ~ 0.46 dex at M* = 10 11 M Q and 
ci og M h ~ 0.7 dex at M* = 10 115 M . As a result, 
the data are particularly sensitive to cri og M, at large 
M*. Therefore, accounting for the mass dependence of 
<?logM» has little impact on the results. 

Finally, we note that our derived values for the param- 
eters of the SHMR should be independent of ci og M» • The 
observables (g-g lensing, clustering, and SMF) do them- 
selves depend on cri og M,but since we account for o-i og M,in 
the model, the extracted SHMR should reflect the true 
underlying physical relationship between halo and stel- 
lar mass (this however would not be the case if one were 
to neglect t7i og M»). In other terms, our model fully ac- 
counts for Eddington bias in all three observables (also 
see discussion in Paper I). 

4.3. Parameters in the model 

To model the central occupation function, we use five 
free parameters (Mi, M*, , /?, 8, 7) to model / SHUH and we 
leave <Ji og M, as an additional free parameter. As de- 
scribed in Paper I, the central occupation function for 
a sample of galaxies more massive than M * is expressed 
as: 



(N ccn {M h \Ml)) = 
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Fig. 4. — The three sources of error that contribute to cr™ gM-B for the SlG_MOD2 model. At all redshifts, o"j™ gM ^ is dominated by model 
error followed by photo-2 error and finally by K-band magnitude error. The rise in the total error at high masses in the low-z bin is 
caused by the model error term. This is likely a result of degeneracies and template errors in the parameter space of models that provide 
reasonable fits to the early-type SEDs that dominate the high-mass end at low redshift. Indeed, the estimated stellar mass PDFs at high 
mass in the low-z bin tend to show more bi-modality. The difference in error across the full mass range from this effect is small however, 
corresponding to only 0.03 dex. 



1 - erf 



logioM - log 10 (/ SHMR (M /l )) 



ClogM, 



(14) 



Our model also has an additional 5 parameters that 
are necessary to model the satellite occupation function, 
(N sa t). For a set of galaxies more massive than threshold 

Mi, (N sat ) is: 



(N sat (M h \MD) -- 
(N ccn (M h \Mt)) 



/sHMR (j^*) 

M sat 



exp 



/shmr(M*) 



(15) 

The free parameters that determine satellite occupa- 
tion as a function of stellar mass are /3 S at, B sat , /3 C ut, 
B cut , and a sat . The first two parameters, /3 sat and B sat , 
determine the amplitude of (N sat ). The second two pa- 
rameters, /3 cu t and B cut , set the scale of the exponential 
cut-off. These parameters enter into (-/V sat ) as: 



10 12 M, 
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and 
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B s 







/sHMR {M* 1 ) 

10 12 M o 



/shmr(^) 



10 12 M, 







Pea. 



A=ut 



(16) 



(17) 



Finally, a sa t represents the power-law slope of the satel- 
lite mean occupation function. We set a sat = 1 for 
all samples which should be a good choice because 
the theoretical expectation is that the number of sub- 
halos above a given mass scales linearly with halo mass 
dKravtsov et all 120041: iConrov et all 120061: iMoster et ail 
120101: iTinker et all \20W . Results from group cata- 
logs also indicate that a sat ~ 1 (jCollister fc Lahavl 
[20051: lYang et all [2009h . Previous HOD analyses of 
clustering results at v arying redshifts also vary little 
from a value of unity (iZehayi et al.l 120051 iZheng et all 
20071: Ivan den Bosch et all 120071: ITinker et all l2007t 



Zehavi et alll2010() . 



In total, our model contains ten free parameters and 
one fixed parameter (a sa t). A summary and description 
of these parameters can be found in Table 5] and also in 
Paper I. 

4.4. Covariance matrices 

The COSMOS survey covers a relatively small volume 
and therefore sample variance effects must be taken into 
account. The volumes probed in the three redshift bins 
are given in Table [1] and vary from 0.88 x 10 6 Mpc 3 for 
the ACS region in z-y (0.22 < z < 0.48) to 4.39 x 10 6 
Mpc 3 for the Subaru region in z a (0.74 < z < 1.0). The 
volumes sampled by COSMOS are too small to obtain an 
accurate estimate of the sample variance from the data 
itself. 

A series of mock catalogs are used to calculate the co- 
variance matrices for all three obscrvablcs. Details re- 
garding the construction of these mocks can be found in 
Paper I. Briefly, COSMOS-like mocks are created from 
a single simulation (named "Consuelo") 420 hr 1 Mpc 
on a side, resolved with 1400 3 particles, and a particle 
mass of 1.87xl0 9 h- 1 MgB This simulation can ro- 
bustly resolve halos with masses above ~ 10 11 /i -1 M 
and is part of the las Damas suitifl (McBride et al. in 
prep). We create mocks for the three redshift intervals: 
2i = [0.22,0.48], z 2 = [0.48,0.74], and 23 = [0.74, 1]. For 
each redshift interval, we construct a series of mocks cre- 
ated from random lines of sight through the simulation 
volume that have the same area as COSMOS and the 
same comoving length for the given redshift slice. This 
yields 405 independent mocks for the z\ bin, 172 mocks 
for the Z2 bin, and 109 mocks for the Z3 bin. For each 
redshift bin, mocks are created from the simulation out- 
put at the median redshift of the bin (see Table [1} . 

We converge on the mocks used for estimating errors 
using an iterative method. To begin with, we find an 
initial best-fit model to the data, without using any co- 
variance matrices. This initial fit is used to populate the 

5 In this paragraph, numbers are quoted for Ho = 100 h 
km s _1 Mpc -1 

6 Further details r egarding this simula tion can be found at 
http : //lss .phy . vanderbilt . edu/lasdamas/simulat ions .html 
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TABLE 4 

Parameters in model 



Parameter 


Unit 


Description 


/AT \ n _ / AT \ 

(iV CBn ) or (-/Vsat) 


T " In ] 

.brec/rixed 


Ml 


M Q 


Characteristic halo mass in the SHMR 


(N ccn ) 


free 


Af*,o 


M 


Characteristic stellar mass in the SHMR 


(N cen ) 


free 


P 


none 


Faint end slope in the SHMR 


(iVcen) 


free 


S 


none 


Controls massive end slope in the SHMR 


(-Ween) 


free 


7 


none 


Controls the transition regime in the SHMR 


(Ween) 


free 


°logM» 


dex 


Log-normal scatter in stellar mass at fixed halo mass 


{N ccn ) 


free 


Asat 


none 


Slope of the scaling of M sa t 


(AT sa t> 


free 


-Bsat 


none 


Normalization of the scaling of M sa t 


(iVsat) 


free 


/Scut 


none 


Slope of the scaling of M cu t 


(iVsat) 


free 


Bent 


none 


Normalization of the scaling of M cut 


(N aat ) 


free 


Qsat 


none 


Power-law slope of the satellite occupation function 


(iVsat) 


fixed 



mocks and to create a first set of covariance matrices. 
We then re-fit the data using these co-variance matri- 
ces and use the best-fit HOD models to create our final 
covariance matrices. All results presented here use this 
final set of covariance matrices. 

5. RESULTS 

We now present the results of fitting the model de- 
scribed in Paper I to the observed COSMOS g-g lensing, 
clustering, and stellar mass functions. 

5.1. Constraining the model 
In order to fit the model to the data, we minimize: 

(18) 

i 3 

where the sum over i and j indicates summation over the 
different stellar mass bins and thresholds, respectively. 
For the SMF and each stellar mass sample in w(9) and 
AS, x 2 lS calculated by: 

X 2 = S(Sn-l/n)C^ 1 (x,- tt ), (19) 
n.l 

where x n is the model calculation of the quantity for data 
point n, y n is the measurement, and C _1 is the inversion 
of the covariance matrix. 

To obtain the posterior probability distributions for 
the parameter set (Mi, M* >0 , /?, S, 7, cri ogM », Asat, -B sa t, 
/?cut) B cn t) we implement a Markov Chain Monte Carlo 
(MCMC) algorithm as follows: 

1. We sample the region of interest around the best 
fit by initializing all chains close to the minimum 
X 2 solution with varying initial random spreads. 

2. We apply conservative limits on all ten parameters. 
As shown in Figure [5J the recovered posterior dis- 
tributions are independent of these limits. 

3. The covariance matrix is updated for the first 3000 
elements in the chain, and then held fixed for the 
remainder of the chain. 

4. We run 6 chains per case, with ~ 40,000 steps. 
We discard all elements before the covariance ma- 
trix is constant (the first 3000), leaving ~ 37,000 
elements per chain. 



5. We use the GetDist package provided with Cos- 
moMC (|Lewis fc Bridle! f2002t) for computing con- 
vergence diagnostics. The chains are tested for con- 
vergenc e and mixing with the Gelman-Rubin cr i- 
tcrion (|Gelman fc Rubinl U99l IGilks et alJll996lh 
We impose a limit on the worst R — 1 statistic of 
R — 1 < 0.01 and observe that for all cases under 
study the mean and marginalized posterior distri- 
butions are well constrained. 

6. We use GetDist output statistics to estimate con- 
fidence limits on all parameters. 

5.2. Fits to the data 

Figures [5l |H and [7] show the best fit models for 
each of the three redshift bins for the galaxy cluster- 
ing, the g-g lensing, and the SMF. The upper panels 
show the angular correlation function w(0) for stellar 
mass threshold samples. The middle right panel shows 
the COSMOS SMF which is measured for all galaxies 
with M* > 10 8 - 7 Mq for zi, M* > 10 9 3 M Q for z 2 , and 
M* > 10 9 - 8 M Q for z 3 . The dotted blue line shows the 
SMF of satellite galaxies from our model. The lower pan- 
els show the g-g lensing signals for the stellar mass bins 
defined in Table [3l When looking at Figures and 
one must keep in mind that the data points are cor- 
related (see Paper I for the covariance matrices) and so 
it is difficult to evaluate "by eye" whether or not the fits 
are adequate. The reduced \ 2 f° r tire fits are 1.7, 1.6, 
and 1.9, for z%, z%, and Z3 respectively. 

For our z\ sample, we have compared the COSMOS 
mass function wi th previously published mass func- 
tions from SPSS jLi fc W hitel [20091: [Baldrv et all [20081: 
iPanter et al.l 120071 ). Because of the use of photometric 
redshifts in COSMOS, we expect a larger stellar mass 
measurement error which should cause Eddington bias 
and lead to an inflated observed SMF at the high mass 
end in COSMOS compared to SDSS. A more in-depth 
comparison of the mass functions is discussed further on 
in §[5H 

5.3. Parameter constraints and redshift evolution 

Table [5] gives the best fit values from GetDist for all 
ten parameters and for the three redshift bins. Figure 
[5] shows the one dimensional and two dimensional joint 
marginalized constraints on the parameters for the z<i 
bin. All of the parameters have well behaved, uni-modal 
distributions. In the interest of brevity, we have not in- 
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10.8<log(M„) b) _ r 
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Log 10 (M» [My) 



10.89<log(M*)<11.12i) 




10.64<log(M,)<10.89j) 




10.3<log(M„)<10.64 k) 




9.82<log(M*)<10.3 1) 



9.2<log(M*)<9.82 m) 



8.7<log(M, ( )<9.2 n) 




1 1(T 0.1 1 10" 

Physical transverse distance R [ Mpc ] 

Fig. 5. — Best fit model for the z\ (0.22 < z < 0.48) redshift bin (blue solid line). Panels a) to f): amplitude of the angular correlation 
function w as a function of angular separation 9 (in arc-seconds) in stellar mass thresholds. Note that in this redshift bin, the amplitude 
of w at large separations is artificially deflated by integral constraint (but this is accounted for in the fitted model). Panel g): COSMOS 
SMF for M« > 10 87 M r (completeness limit for this redshift bin). For ref e rence , in panel f), we also show the S PSS mass fu n ctions 
from ILi fc White! (120091 1 (triple schecter fit, black dashed line), iBaldrv et all 1120081 ) (green dash-dot line) and from IPanter et all (1200711 
(orange, dash-dot-dot line). The dotted blue line in panel g) shows the SMF of satellite galaxies for the best fit model. Panels h) to n): 
galaxy-galaxy lensing signal in stellar mass bins. Note that in panel m), there is a negative data point represented by a grey square (this 
can occur due to noise and is not a concern). The lensing signal is decomposed into four components: the baryonic term (red dotted line), 
the central one-halo term (green dashed line), the satellite one-halo term (orange dash-dot), and the two-halo term (grey dash-dot-dot-dot 
line). 
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Fig. 6. — Same as Figurc[5]but for the second redshift bin. Z2- Panels a) to e): angular correlation function. The angular correlation 
function is less affected by integral constraint in this redshift bin since the volume probed is 2.3 times larger than Z\. Panel f): COSMOS 
SMF for M* > 10 9 ' 3 Mq. The dotted blue line in panel f) shows the SMF of satellite galaxies for the best fit model. Panels g) to m): 
galaxy-galaxy lensing signal. 
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Fig. 7. — Same as Figure \E\ but for the third rcdshift bin, 23. Panels a) to d): angular correlation function. The angular correlation 
function is less affected by integral constraint in this redshift bin since the volume probed is 5 times larger than z%. Panel e): COSMOS 
SMF for M* > 10 9 8 M . The dotted blue line in panel e) shows the SMF of satellite galaxies for the best fit model. Panels f) to k): 
galaxy-galaxy lensing signal. 



The stellar-to-dark matter connection in COSMOS 



15 



TABLE 5 

Best fit parameters for the 3 redshift bins 



Parameter 


2i = [0.22,0.48] 


z 2 = [0.48, 0.74] 


z 3 = [0.74, 1] 


SIG.MOD1: 


log 10 (Mi) 


12.520 ±0.037 


12.725 ±0.032 


12.722 ±0.027 


log 10 (Af*,o) 

/3 


10.916 ±0.020 


11.038 ±0.019 


11.100 ±0.018 


0.457 ±0.009 


0.466 ± 0.009 


0.470 ± 0.008 


5 


0.566 ±0.086 


0.61 ± 0.13 


0.393 ± 0.088 


7 


1.53 ±0.18 


1.95 ±0.25 


2.51 ±0.25 


°"logM„ 


0.206 ±0.031 


0.249 ±0.019 


0.227 ±0.020 


B C ut 


1.47 ±0.73 


1.65 ±0.65 


2.46 ± 0.53 


Bsat 


10.62 ±0.87 


9.04 ±0.81 


8.72 ± 0.53 


/Scut 


-0.13 ±0.28 


0.59 ±0.28 


0.57 ±0.20 


Asat 


0.859 ± 0.038 


0.740 ±0.059 


0.863 ±0.053 



SIG.MOD2: 



log 10 (Mi) 

lo Sio( M *,o) 

/8 
S 

7 

i * 
°"logM„ 
Bcut 
Bsat 
/3cut 
Asat 



12.518 ±0.038 
10.917 ±0.020 
0.456 ± 0.009 
0.582 ±0.083 
1.48 ± 0.17 
0.192 ±0.031 
1.52 ± 0.79 
10.69 ±0.89 
-0.11 ±0.29 
0.860 ± 0.039 



12.724 ± 0.033 
11.038 ±0.019 
0.466 ± 0.009 
0.62 ± 0.12 
1.93 ± 0.25 
0.245 ± 0.019 
1.69 ± 0.65 
9.01 ±0.81 
0.60 ± 0.27 
0.740 ±0.059 



12.726 ±0.028 
11.100 ±0.017 
0.470 ± 0.008 
0.47 ±0.10 
2.38 ± 0.24 
0.220 ±0.019 
2.57 ±0.56 
8.66 ±0.53 
0.58 ± 0.20 
0.863 ±0.053 



In the SIG.MOD2 case we fit for a 



loglvt 



whereas in the SIG.MOD1 



case we fit for cri og M* 

Note: the quoted errors are always marginalized over all other 
parameters. 

eluded equivalent Figures for z\ and Z3 but they are sim- 
ilar to Figure [51 All parameters are well constrained in 
the three redshift bins. 

Table [5] lists the best fit values for SiG_MODl where we 
have assumed that cri og M, is constant and for SIG_mod2 
where we have explicitly accounted for stellar mass de- 
pendent errors. We find no strong difference in our re- 
sults, regardless of which model we adopt. We conclude 
that our results are robust to the effects of mass depen- 



dent scatter. In the SIG_mod2 case, we model a 



logM, 



and 



assume that CTi og M» is the sum in quadrature of tTj ogMt 
(which is assumed to be constant) and cr^gM, • Therefore, 
the quantity that we actually fit for in the SIG_MOd2 case 
is 0j ogMit (whereas in SlG_MODl we fit for ai og M, )■ This is 
why, as expected, the best-fit scatter in Table[5]is slightly 
lower for SlG_MOD2 compared to SIG_mod1. Note that 
we are not claiming to actually extract meaningful values 
for the intrinsic scatter in stellar mass at fixed halo mass. 
To do so would require a more thorough error analysis 
which is beyond the scope of this paper. Overall, our 
conclusions regarding (Ti og M» are twofold. Firstly, we can 
safely assume that cri og M» is constant and ignore any mass 
dependent effects induced by measurement error. This 
is due to the fact that all three observables are primarily 
sensitive to the effects of <7i og M„ a t large M* where the 
slope of the SHMR increases sharply. Similar conclusions 
were reached by BIO who find that the effects of scatter 
are insignificant below M* = 10 10 5 M Q where the slope 
of the SHMR is not steep enough to have a significant 
impact. Secondly, we find that (7\ og M r ~ 0.23 ± 0.03 dcx, 
in broad agreement with previous results. For example, 



B10 estimate that c™ gM<i =0.07 dex and crf ogM> =0.16 
dex, resulting in a total scatter of <7i og M, = 0.175 dex 
(their total scatter is lower than ours as expected be- 
cause we have a larger o"i™ gM component). 

Figure [9] shows the measured redshift evolution for all 
ten parameters. Previous work on this topic has been 
limited by systematic differences in stellar mass estimates 
between low and high-z surveys which can be of order 
0.25 dex (excluding IMF uncertainties) according to BIO. 
We stress that the results in this paper have been derived 
in a homogeneous fashion from high to low redshift. Our 
conclusions regarding redshift evolution should thus be 
robust in this respect. 

The most striking result from Figure [HI is the redshift 
evolution in the two parameters M\ and M*^. This is 
one of the major results of this paper which we will dis- 
cuss in more detail in the following section. Apart from 
these two parameters, there is marginal evidence for some 
evolution in 7 and B sat - No strong evolution is detected 
in the remaining six parameters. It is interesting to note 
that j3 (which controls the low mass slope of the SHMR) 
remains constant at j3 ~ 0.46. We will provide an inter- 
pretation of this result in the discussion section. 

5.4. The SHMR from z = to z = 1 

Figure [T0l shows the best fit SHMR for Z\ compared to 
a variety of low red s hift con st raints from weak l ensing 
(iMandelbaum et all l2006bt iLeauthaud et all 120101: 



Hoekstr al 120071). abun d ance matching (jMoster et al.1 



2010t Behroozi et al.l I2010D. sate llite kinematics 
i Conrov et a ll 120071: iMore et all f2010h and the Tully- 
Fisher relation dGeha et al.l [2006T iPizagnol [20061: 
ISpringob et "all 120051 : iBlanton et al.1 120081 ) (see § [H 



for a more in-depth comparison). Most of the data 
are in broad agreement and show clear evidence for 
a variation in the dark-to-stellar mass ratio with a 
minimum of M h /M t ~ 27 at M* ~ 4.5 x 10 10 M and 
M h ~ 1.2 x 10 12 M Q . As demonstrated by B10, however, 
meaningful and detailed comparisons between various 
data-sets are hampered by systematic uncertainties in 
stellar mass estimates. For this reason, we will mainly 
focus in what follows on the conclusions that can be 
drawn by inter-comparing the three COSMOS redshift 
bins. 

At low stellar mass, Mh scales as Mh oc (M,) 046 and 
this scaling does not evolve significantly with redshift 
from z = 0.2 to z = 1. At high stellar mass, the SHMR 
rises sharply at M* > 10 10 5 M Q as a result of which 
ciogM h (the scatter in halo mass at fixed stellar mass) also 
increases and M* is clearly no longer a good tracer of halo 
mass. For example, a galaxy with ~ 10 11 " 3 Mq may 
be the central galaxy of group with Mh ~ 10 13 — 10 14 M Q 
or may also be the central galaxy of a cluster with Mh > 
10 15 M@. 

A quantity that is of particular interest is the mass at 
which the ratio Mh/M* reaches a minimum. This min- 
imum is noteworthy for models of galaxy formation be- 
cause it marks the mass at which the accumulated stellar 
growth of the central galaxy has been the most efficient. 
We describe the SHMR at this minimum in terms of the 
"pivot stellar mass", Mf lv , the "pivot halo mass", M^ lv , 

and the "pivot ratio", (M h /M*)P iv . Note that M? iy and 

M£ 1V are not simply equal to Mi and M*. . Indeed, the 
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Fig. 8. — One and two dimensional joint-mean and marginalized distributions for all ten parameters for the Z2 bin. Dotted lines in the 
ID distributions (and shaded contours in the 2D planes) show mean likelihoods of samples and solid lines show marginalized probabilities. 
In the 2D planes, blue to red denotes the region of lowest to highest likelihood density. Solid contours in the 2D plane represent the 68% 
(lcr) and 95% (2cr) confidence regions. 



mathematical formulation of the SHMR is such the pivot 
masses depends on all five parameters. The three param- 
eters that have the strongest effect on the pivot masses 
are Mi, M* q, and 7 (see Paper I). 

Figure [11] shows the redshift evolution of the SHMR. 
Three points are of particular interest in this Figure. 
First, we detect no strong redshift evolution in the low 
mass slope of the SHMR (M» < 10 10 2 M ). Indeed, as 
highlighted in the previous section already, (3 is remark- 
ably constant out to z = 1. In Paper I we have shown 
that f3 regulates the faint end slope of the SMF so in 
other terms, we could also state that the faint end slope 
of the SMF shows remarkably little redshift evolution. 



We do however find that the amplitude of the low-mass 
(M* < 10 10 - 2 M ) SHMR was higher at earlier times so 
that galaxies at fixed stellar mass live in more massive 
halos at earlier epochs. We will return to this result in 
the discussion section. 

Second, we detect a strong redshift evolution in M* 1V 
and M^ IV which is shown more explicitly in Figure 1121 
The detected evolution is such that both the halo mass 
and the stellar mass for which the accumulated stellar 
growth of the central galaxy has been the most efficient is 
smaller at late times than at earlier times. This trend is a 
manifest ation of at least one meaning of the term "down- 
sizing" (jCowie et all [19961 iBrinchmann fc Ellisl 120001 : 
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Fig. 9. — Redshift evolution in the ten parameters. Red diamonds indicate evolution in the five parameters that describe the SHMR 
(Mi, M t fi, /3, 8, 7). Orange squares show the scatter in stellar mass at fixed halo mass (o"iogM* )• Blue triangles show the remaining four 
parameters that describe the satellite occupation function (/3 sa t, B Ba ,t, /3cut, B cu t). Grey data points shows the best fit parameters when 
we assume a stellar mass dependent scatter due to photometric redshifts (the SIG-MOD2 model) . 



I Juneau et al.l 12005). Originally, the term downsizing re- 
ferred to the notion that the maximum K-band luminos- 
ity of galaxies above a sp ecific star formatio n rate thresh- 
old decreases with time (jCowie eFaDEHl). Since then, 
downsizing has been widely employed to more generally 
describe the behaviour in which a certain parameter that 
regulates galaxy formation decreas es with time (for re- 
cent discussions on downsizin g see iFontanot et all 120091 
and lConrov fc Wech slcr 200^). Our results show strong 
evidence for downsizing in both the pivot halo mass and 
the pivot stellar mass. We have already remarked in 
the previous section that a strong evolution in Mi and 
M* is seen in Figurc[9] Although these two parameters 

are not strictly equal to M|" v and M^ IV , they do have a 
strong impact on the location of the pivot masses. Thus, 
the evolution seen in Mi and M*^ is directly related to 
the observed downsizing behaviour in the pivot masses 
that is apparent in Figures [11] and [T^l The pivot stel- 
lar mass evolves from M* piv = 5.75 ± 0.13 x 10 10 M Q at 
z = 0.88 to Mf iv = 3.55 ± 0.17 x 10 10 M Q at z = 0.37 
with an evolution detected at lOc. We note that all er- 
rors have been derived by marginalizing over all other 
parameters. 

The evolution in M|" v varies smoothly in the three red- 
shift bins, however the evolution in M plv is less smooth, 
in particular in the Z2 bin. We suggest that M£ 1V is more 

sensitive to sample variance than M* IV . Indeed, the first 
order effect of sample variance is to change the normal- 
ization of the SMF (see Paper I) ; this will directly affect 



Mi and thus Ml 



In summary: M pl 



is sensitive to 



sample variance between redshift bins whereas M* IV is 
sensitive to systematic errors in stellar mass measure- 
ments between redshift bins. 

Finally, at high masses (M„ > 10 11 M Q ) there is an 
interesting hint that the amplitude of the SHMR is de- 
creasing at higher redshifts, but we lack the statistics for 
a clear detection, mainly due to the small volume probed 
by COSMOS. 



5.5. Comparison with previous work 

Figure [TU] compares our z\ SHMR to previous work 
on this topic at low redshift. The general picture that 
emerges from Figure[l0]is one of remarkable broad agree- 
ment between various methods on the overall shape of 
the SHMR. In detail however, meaningful comparison 
between various surveys are severely limited by system- 
atic differences in stellar mass estimates (^0.25 dex be- 
tween different surveys). For this reason, we mainly focus 
on qualitative comparisons in this section. All results 
have been adjusted to our assumed value of Ho = 72 
km s _1 Mpc -1 and unless stated otherwise, halo masses 
are converted to Myn o h assu ming a NFW profile and a 
iMufioz-Cuartas et al.l (|20100 mass-concentration relation 
for a WMAP5 cosmology when necessary. All results 
quoted here assume either a Kroupa or a Chabrier IMF. 
Since the systematic shift in M* between these two IMFs 
is small (~ 0.05 dex), we do not adjust for this difference. 
We also do not attempt to correct for differences in the 
assumed cosmological model. 

5.5.1. Comparison with previous work: low redshift 

iMandelbaum et akl (|2006al ) have used g-g lensing in 
the SDSS to measure halo masses for len s galaxies at 
z ~ 0.1. Since Mandel baum et akl (|2006al) only present 
their results as a function of both M* and color, the data 
points in Figure [TOl have been re-computed as a function 
on M* only (Rachel Mandelbaum, priv comm). Except 
perhaps for one data point at low M*, these results are 
in good agreement with ours. 

At high masses, an alternative method to probe the 
central SHMR is to directly compare the halo masses of 
groups and clusters of galaxies to the masses of their 
central galaxies. Since we are primarily interested in 
<J> c (M*|M/j), it is critical, as much as possible, to use 
halo mass selected samples of groups and clusters for this 
type of comparison. Using samples selected on the ba- 
sis of the stellar mass of the central gal axy, for example, 
would result in biased conclusions. In iLeauthaud et al.1 
(|2010|) we presented a sample of X-ray groups (M 2 oofc ~ 
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Fig. 10. — Top panel: Inferred SHMR in the z\ redshift bin compared to other low redshift measurements from weak lensing (WL), 
abundance matching (AM), satellite kinematics (SK), and the Tully Fisher relation (TF). The COSMOS zi results are shown by the solid 
dark blue line a nd the shaded grey region represents the one sigma error on the SHMR. This SHMR represents (log 10 (M* (M h ))). With 
the exception of lConrov et al. (2007), all data points either use or have been converted to this same averaging system . Overall, there is a 



broad agreement between various probes. Detailed comparisons between various data-sets however, are limited by systematic differences in 
stellar mass estimates due to varying assumptions (e.g., star formation histories, extinction laws, stellar population models). Bottom panel: 
Dark-to-stellar mass ratio as a function of stellar mass. We observe a clear variation in A^oot/Af, with A/2006/ 'M* reaching a minimum of 
M h /M* ~ 27 at M* ~ 4.5 X 10 10 M and M 200b ~ 1.2 X 10 12 Mq. The dark-to-stellar mass ratio rises sharply at Af» > 5 X 10 10 M so 
that a cluster of halo mass A^oob ~ 10 15 Mq will reach a ratio of M2oob/A4» ~ 2000. Note that this ratio only refers to the ratio between 
the halo mass and the stellar mass of the central galaxy. For example, in the case of clusters, we are comparing the ratio of the cluster halo 
mass to stellar mass of the central Brightest Cluster Galaxy (BCG). 



10 13 - 10 14 M ) in COSMOS for which we have cal- 
ibrated the relationship between halo mass and X-ray 
luminosity (Lx) using g-g lensing. The expected scat- 
ter in halo mass at fixed Lx is of order 0.1 3 dex so the 
sample presented in lLeauthaud et all (|2010T l is halo mass 
selected to a good approximation. In parallel, George et 
al. (in prep) have constructed an algorithm to identify 
the central galaxies of these groups and have used the 
weak lensing signal itself to optimize the algorithm by 
maximizing the weak lensing signal at small radial sep- 
arations from the central galaxy. The grey squares in 
Figure [TU] report the stellar mass of the central galaxy 
versus M/j for groups at 0.22 < z < 0.48 and with a 
high quality flag. These data points are directly compa- 
rable to ours since we have used exactly the same stellar 



masses and confi r m tha t our results are consistent with 
iLeauthaud et all (|2010l) . 

We present a similar exercise for a sample of X-ray lu- 
minous clusters (A68, A209, A267 A383, A963, A 1689, 
A1763, A2218, A2390, A22 191 fromlHoekstral (120071) with 
weak lensing masses from iMahdavi et all (|2008fT The 
centra l galaxies of these clu sters have been studied in de- 
tail bv lBildfell etaLl (|2008D . Using the same stellar mass 
code and assumptions as in this paper, we have com- 
puted stellar masses for the central cluster galaxies using 
a compilation of optical data provided by Chris Bildfcll. 
The results are shown by the red asterisk points in Fig- 
ure [101 Unfortunately, these mass estimates are based on 
just two optical bands (B-band and R-band) and as such 
will have larger uncertainties than the COSMOS stellar 




Fig. 11. — Left panel: Redshift evolution in the SHMR. Because the COSMOS results have been derived in a homogeneous fashion, we 
can inter-compare the three COSMOS redshift bins. The low mass slope of the SHMR remains constant with ft ~ 0.46. We detect a clear 
evolution in Mf lv and M? 1V which are respectively the stellar and halo mass at which M20Qb/M t reaches a minimum (see right panel). 
At high masses, there is an interesting hint that the amplitude of the SHMR is decreasing at higher redshifts, but we lack the statistics to 
make a clear detection, mainly due to the small volume probed by COSMOS. We have also plotted the SHMR reported by BIO. However, 
caution must be taken when making a direct comparisons between COSMOS and BIO because our stellar masses have been derived under 
different assumptions. According to BIO, the level of systematic uncertainty in stellar masses is of order 0.25 dex. Given this 0.25 dcx 
systematic uncertainty, we are in broad agreement with BIO but a direct comparison would require a more homogeneous analysis between 
COSMOS and SDSS. Flight panel: redshift evolution in Afgoob/-^*- Both M# lv and M^ lv exhibit downsizing trends, decreasing at later 
epochs. This effect is shown more explicitly in Figure [12] 



masses used in this paper which are constrained with 
many more filters and normalized to a near-IR luminos- 
ity. We estimate that an additional 0.25 dex stellar mass 
uncertainty should be included when interpreting these 
data points, which may account for their scatter in stellar 
mass seen in Figure [TU] With this cautionary note, plus 
the additional caveat that this sample is not as homo- 
geneo usly selected as the groups from iLeauthaud et al.l 
(|2010t l. and that COSMOS is too small to probe over- 
densities of these masses, the results are nonetheless in 
good a greement with the ex trapolation of our z\ SHMR. 

Both lMoster et ahl ([2010D and BIO have presented con- 
straints on the SHMR and its redshift evolution by us- 
ing the abundance matching technique. A more detailed 
comparison with their work is presented in § 15.5.31 and § 

IConrov et al.l (|2007t) have used the kinematics of satel- 
lite galaxies in SDSS and DEEP2 (Davis 2003) to probe 
the SHMR at z ~ 0.06 and at z ~ 0.8. Their low redshift 
results are shown by the blue diamonds in Figure [TU] A 
30% downward correction to halo masses has been ap- 
plied due to incompleteness effects as described in their 
paper. In addition to systematic differences in stellar 
mass estimates, direct comparisons with our results are 
further complicated by the fact that our mod el descr ibes 
(log 10 (M*(M /l ))) whereas the IConrov et all (|2007l ) re- 
sults represent (M^(M*)). The two averaging systems 
will yield different results: (M/i(M*)) will be increas- 
ingly biased low with respect to (log 10 (M*(M/j))) with 
increasing a\ og m, and for steeper values of the slope of 
t he SHMR. 

iMore et al.l (|2010l) have used SDSS data to probe the 
halo masses of ~ 3900 central galaxies in the range 



0.02 < z < 0.072 using the kinematics o f satellite 
galaxies. In their paper, IMore et al.l (pOlOt) have an- 
alyzed red and blue galaxies separately. For Figure 
[TU1 we have asked the authors to provide the data for 
all central galaxies function of M*, irrespective 

of color, and also to convert their results to reflect 
the mean log-relation ((log 10 (M*(M/j))) as opposed to 
(Mh(M*))). Overall , there is some d isagreement be- 
tween our results and IMore et al.l (|2010T l regarding of the 
general shape of the SHMR. Indeed, our results display 
a more strongly v arying power-law index compared to 
IMore et al.l (|2010l ). This disagreement is perhaps more 
apparent in the lo wer panel of Figure [KB Indeed, the 
IMore et al.l (|2010D results display a fairly broad mini- 
mum in AhoQb/M* whereas our results predict a more 
strongly varying M2006 /M* ratio that reaches a minimum 
at M* - 4.5 x 10 10 M Q . IMore et al.l (|2010D have sug- 
gested that satellite kinematics may yield halo masses 
that may be systematically higher by a f actor of 2-3 
than o ther methods at low M». Lowering the IMore et al.1 
(|2010[) results at low M* would certainly bring their re- 
sults into better agr eement with ours in terms of the 
shape of the SHMR. IMore et al.l (|2010l) also provide es- 
timates for ciogM,- They find aiogM, = 0.19°:q 3 03 for red 
centrals and cri og M„ = 0.15^g 2 07 for blue centrals. Both 
of these values are in broad agreement with our estimate 
of uiogM, ~ 0.23 dex. 

In the z\ redshift bin, the COSMOS results are lim- 
ited by completeness to Af* > 10 8 7 M Q . Nonetheless, 
it is of interest to see how our results extrapolate to 
galaxies of even lower stellar masses, even though mea- 
surements of Mh for such low mass galaxies are fraught 
with difficulties and for the most part limited to the lo- 
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Fig. 12. — Redshift evolution in the pivot halo mass, the pivot stellar mass, and the pivot ratio. Dark blue circles show the results from 
this paper. We detect a striking evolution in Mj? lv with lv decreasing at later epochs. This downsizing in Mf lv is accompanied by a 
downsizing in M^ lv with M plv also decreasing at later epochs. We detect no strong evolution in the pivot ratio which remains constant at 
(M^/M*) plv ~ 27. Green diamonds: comparison with B10. Yellow dash-dot line: comparison with IMoster et all ( f20Totl (from Table 7 of 
their paper). The downsizing trend in M^ lv (which is not sensitive to systematic errors associated with M„) is detected by all three studies. 
The normalization and evolution of M^ lv is similar for COSMOS and for B10. However, the normalization of M^' v from Moster et al. 
1120101 ) is higher than our prediction. This could be due to differences in the adopted parametric form of the SHMR. The three studies 
show different results for the evolution of M* lv . We stress that stellar mass estimates computed under different assumptions may present 
relative systematic errors of order 0.25 dex. This error will affect quantitative comparisons for the normalization M* lv (and by consequence, 
also (M h /M*) plv ). Since our results have been derived homogeneously from z = 0.2 to z = 1, our conclusions regarding the evolutionary 
trends of A/f lv should be robust. In contrast, B10 and [Moster et al. (20l7j) infer evolutionary trends from a combination of SDSS data 

and other distinct surveys at higher rcdshifts. Therefore, their results have large systematic uncertainties with res pect to Mf". These 
systematic errors are reflected in the B10 error bars (B10 account for systematic errors and sample variance). The IMoster et al. 1120101) 
result s do not account for systematic error or sample variance from mock catalogs. Adding these two sources of error to the IMosTeretHd] 
1120101 ) results would lead to errors on the pivot quantities of roughly the same order as B10. Given these considerations: both the B10 and 

the Moster et al. (2010) results (provided larger errors are included) are consistent with our detection of an evolving Af* 1 and a constant 
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Fig. 13. — An independent validation of the evolution of the pivot masses using the abundance matching method of B10. Blue circles 
show our detected evolution of the pivot quantities. Magenta diamonds show the evolution of the pivot quantities when the abundance 
matching method of B10 is applied to the COSMOS SMFs. We find that the analysis of B10 applied to the COSMOS results fully agrees 
with our claims concerning the evolution of the pivot quantities. We conclude that the detected downsizing behaviour of M? lv and Mj lv 
is robust to the methodology that is employed. 



cal volume. iBlanton et al.l (|2008|) have presented an ef- 
fort to address the very low mass SHMR (see their Fig- 
ure 12) using measurements of the maximum circular 
velocities from Hi disks around iso lated nearby dwarf 
galaxies. Since IBlanton et al . (2008) have applied crite- 
ria to specifically select isolated galaxies, their sample 
should be dominated by central galaxies and so com- 
parable to our Figure 1101 Upon request, the authors 
provided us with th e full data set from Figure 12 in 
IBlanton et all (|2008D which we have rep roduced in Fig- 
ure [Tpl inc luding addition al data from ISpringob et al.l 
(|2005D and IPizagno et"all (|2007t) based on H T and Ha 



rotation curves respectively. This data compilation is re- 
stricted to galaxies that are isolated and with axis ratios 
b/a < 0.5 in order to minimize inclinatio n uncertainties 
and ex tinction corrections. Halo masses in IBlanton et al.l 
(|2008f ) have been estimated by assuming that the optical 
circular velocity, V op t, is equal to V maX j the maximum 
circular velocity for a NFW halo. V mllx is then con- 
verted to t he viral velocity, V20 0, using N-body calibra- 
tions from iBullock et ail pOOlf ). For galaxy mass halos, 
V max /V2oo — 1-1 — 1-2 under the assumption of no adia- 
batic contraction of the dark matter due to galaxy for- 
mation. When incorporating adiabatic contraction into 
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Tully-Fisher analyses, IGnedin et al.l (|2007t) find a factor 
of ~ 2.5 decrease in the inferred halo mass at fixed stel- 
lar mass. Such a correction would put the Tully-Fisher 
constraints into better agreement with our results. 

To first order, there is relatively go od agreement 
betwe en our SHMR and the data from IBlanton et alJ 
(|2008|) . albeit with a much larger scatter in the Tully- 
Fisher based SHMR than predicted by our results. 
In particular, the dwarf galaxy data points from the 
IGeha et alJ (|2006[ ) sample are in good agreement with 
the extrapolation of our SHMR to lower masses. At 
10 9 M Q < M* < 10 11 M , however, ther e may be some 
indica tion that the halo-masses inferred bv lBlanton et ahl 
(|2008l) arc too high on average compared to our results 
with the possible implication that the V op t/V2oo is larger 
than 1.1 — 1.2. The T4pt/V2oo ratio contains informa- 
tion about the relative importance of baryons versus dark 
matter on galaxy scales: V opt /V2oo > 1.1 — 1.2 would im- 
ply that the baryons have modified the dark matter pro - 
file in the very inner halo regions. IDutton et all ([20lof ) 
have used the IPizagno et al.l (|2007[ ) data in combination 
with a compilation of prior work on the SHMR to place 
constraints on Km t/V200 for late type galaxies (the re- 
quirement that the lPizagno et all200 71 galaxies have suf- 
ficiently extended Ha emission to yield a useful rotation 
curve implies that this is primarily a late type sample). 
They find that V^ p t ~ V20O) however, the normalization 
of the SHMR tha t they employ is more similar to the 
IMore et aLl (|2010D results than t o our s. Therefore, a 
sim ilar analysis as Dutton et al.1 (|2010t) but applied to 
the IPizagno et al.1 (|2007f ) data in combination with our 
results would yield a higher V^pt/V^oo- 

The agreement between th e extrapola t ion o f our 
SHMR to lower masses and t he IGeha et all (12001 sam- 
ple is encouraging. In fact, iBusha et al.1 (j201Q[ ) found 
that a similar scaling continues to work down to the 
faintest satellites of the Milky Way. However, Figure 
[TU1 clearly reveals a lack of data at M* < 10 9 M due 
to the stellar mass completeness limits of current op- 
tical and IR surveys. Pushing the SHMR down to 
10 8 M© < M* < 10 9 M Q using techniques such as de- 
scribed in this paper is clearly an exciting avenue to ex- 
plore and will be facilitated by upcoming very deep op- 
tical and IR surveys such at UltraVista and the Hyper 
Suprime Cam (HSC) survey on the Subaru telescop^ll- 

5.5.2. Comparison with previous work: high redshift 

IHevmans et al.1 (|2006aD have used g-g lensing to esti- 
mate halo masses for a sample of 626 galaxies with M* > 
10 10 - 5 M Q and with 0.2 < z < 0.8 f ro m the 0.25 deg 2 
HST/GEMS survey (|Rix et all 120041 ). IHevmans et~atl 
(|2006a| ) find M v j r /M* = 53+*g at a mean stellar mass 
of M» = 7.2 x 10 10 M . Converting their result to our 
assumed value of H and to M 2 oo& yields M 20 o&/M» ~ 
58^Jy. In a similar mass and redshift range, our re- 
sults produce M200&/M* ~ 34. We note, however, that a 
careful comparison between our work and lHevmans et al.1 
(|2006a[) is limited by several differences in the way the 

7 For example, the UltraVISTA survey of the COSMOS field will 
obtain IR imaging to Y=26.7, J=26.6, H=26.1, K s =25.6, pushing 
low redshift stellar mass completeness limits to below M* = 10 
Mq at z < 0.3. The HSC intermediate layer survey will cover 20 
deg 2 to g=28.6, r=28.1, i=27.7, z=27.1, and Y=26.6 



analyse s have been performed. Firstly, IHevman s et al.1 
()2006al) fit a NFW profile to the g-g lensing signal and 
so the masses that they measure will reflect the mean 
halo mass at fixed stellar mass, which is different than 
our averaging system. Secondly, IHevmans et al.1 (|2006a[ ) 
do not account for the contribution of satellite galaxies 
to the g-g lensing signal and so they will tend to overes- 
timate halo masses. The lens sample of IHevman s et all 
(2006a) is roughly similar to our g-g bin4 in the Z2 red- 
shift range: the contribution to the g-g lensing signal 
from satellites for this sample is shown in panel j of Fig- 
ure |6l Given these two caveats our results are in fairly 
good agreement. 

5.5.3. Comparison with Moster et al. 

IMoster et al.l (|2010( ) derive constraints on the redshift 
evolution of th e SHMR by a bunda nce matching to the 
SDSS SMF of IPanter et al.l ((20071 ) at low redshift and 
to mass functions from the MUNICS survey (detection 
limit K~19.5, area 0.28 deg 2 ) (jDrorv et al.ll2004l ) and the 
GOODS-MUSIC sample (detection l imit K~23.5, area 
143.2 arcmin 2 ) (jFontana et al.1 120061 ) at high redshift. 
There are two main differences between IMoster et al.1 
(l20T0l) and our work. Firstly, we adopt the functional 
form advocated by B10 for the SHMR which is sub- 
expon ential at high M*. In contrast, the IMoster et al.l 
(|2010( ) parametrization asymptotes to a power-law at 
high M*. According to B10, such a parametrization may 
be problematic. Indeed, because the logarithmic slope 
of the SHMR increases with increasing M*, the best- fit 
power-law for high mass galaxies will depend on the up- 
per limit in the avai lable data for the SMF. Secondly, the 
IMoster et all (pOlO ) errors do not reflect po ssible system- 
atic er r ors in stellar mass e stima tes betweenlPanter et al.1 
([2007T ). lDrorv et al.l (|200l . and lFontana et al.1 (|2006f K 

Although we find the same qualitative behaviour as 
IMoster et all (I2010D : M h /M* decreases to a minimum at 
Mh ~ 10 12 Mq and then rises at higher masses, interest- 
ingly, our results differ regarding the evolutionary trends 
of the SHMR. The two parameters for which our conclu- 
sions differ in particular are M* lv and (Mh/M*) piv . In 
Figure [T2] (yellow dash-dot line) we show the evolution 

of M piv , Mf , and (Afj M,)^, as inferred from Table 



7 in IMoster etaJ] (|201C 
Our results agree with IMoster et al.1 (|2010f) in terms of 

the qualitative downsizing trend seen for M piv . However, 
it is interesting to note that our measurements differ with 
respect to the normalization of -M piv . The exact origin of 
this discrepancy remains unclear. In light of the results of 
B10, we hypothesize that this discrepancy may be caused 
by the difference in the assumed parametric form of the 
SHMR. In any case, further investigation regarding the 
source of this discrepancy, though beyond the scope of 
this paper, is clearly warranted. 

Our conclusions differ with respect to M piv and 
(M/,/M*) piv . Whereas our results suggest that M? iv in- 
creases with redshift an d that (M/,,/My) piv remains con- 
stant, in contrast, the IMoster et al.l (|2010D results im- 
ply that M PIV is constant with redshift and that instead, 
(Mh/M*) plv increases with redshift. We hypothesize that 
th is discrepancy is sim ply due to the fact that the errors 
in IMoster et al.l l)2010h are likely to be under-estimated. 
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Indeed, accounting for sample variance with mocks as 
well as for syst ematic differences i n the relative stella r 
masses betweenlPanter et al.l p007l) . iDrorv et all (|200l ). 
iFont ana eTaTl (|2006ft would lead to similar errors as BIO 
(Figure fT2l green diamonds). Indeed, the evolution that 

we detect in M* lv is ~0.21 dex from z = 0.37 to z = 0.88 
which is similar to the expected error in stellar mass es- 
timates between different surveys. This could perhaps 
also explain why we reach similar conclusions regarding 
the evolution of Mh (which should be less affected by 
systematic errors associated with M*) but not M* lv or 

{M h /M*y™. 

5.5.4. Comparison with Behroozi et al. 

The closest comparison with our work is BIO since 
we employ the same functiona l form for the SHMR , the 
same halo mass function from iTinker et all (|200l) . and 
we both account for the effect of scatter in the SHMR 
and for sample variance in the data using mock catalogs. 

BIO derive constraints on the redshift evolution of the 
SHMR by a bund ance matching to the SDSS SMF of 
iLi fe White! (|2009) at low redshift and to mass functions 
from the FIDEL Legac y Project in the extended G roth 
strip at higher redshifts flPerez-Gonzalez et al.ll2008h . As 
a result of the fact that BIO combine data from distinct 
surveys, their systematic uncertainties on the evolution 
of the SHMR are fairly large. We also note that there 
are differences between the stellar mass estimates used 
in BIO and in this paper (see § I5.6[) which will lead to 
normalization differences in Af* IV and (M/ l /M„») plv for 
example. 

Another difference between BIO and our work is the 
treatment of satellite galaxies. In our model, the SHMR 
only applies to central galaxies and satellites are mod- 
elled via (-/V sa t). Indeed, we require a more sophisticated 
treatment of satellites in order to fit the clustering and 
the g-g lensing for which the satellite term plays a larger 
role than in the SMF (the satellite term is sub-dominant 
at all scales for the SMF). In contrast, BIO assume that 
the SHMR applies also to satellite galaxies, on condi- 
tion that the "halo mass" for satellite galaxies is defined 
as the halo mass at the epoch when satellites were ac- 
creted onto their parent halos (the "infall mass" , M in f a n). 
Thus, there could be subtle differences between the two 
methods due to the treat ment of satellite gal axies (for 
example, see discussion in iNeisteTn et al. 201l|). 

Figure [T^] (green diamonds) shows the prediction from 
BIO for the pivot quantities. Our results are in strik- 
ing agreement with BIO with respect to the evolution 
of M£ 1V . The errors from BIO are larger for M* IV 
and (M h /M r ) piv . Thus our results agree with BIO in 
terms of qualitative evolutionary trends for Mf lv and 
(Mh/M*) plv . There is a normalization difference for 
M* IV between BIO and our results. However, this nor- 
malization offset is not unexpected given systematic dif- 
ferences due to varying assumptions for stellar mass es- 
timates. We address this issue further below. 

5.5.5. Testing the redshift evolution of the pivot masses 

Given the striking downsizing signal that we detect for 
the pivot masses, we would like to check that this result 
is not an artifact of the method we are using. As such 



we have also re-analysed the three COSMOS stellar mass 
functions using the abundance matching method of BIO. 
In this test, the method followed here is identical to that 
of BIO but with different constraints on the stellar mass 
function. Specifically, we use the three cosmos SMFs to 
constrain the redshift evolution of the SHMR and do not 
include any SDSS data. The results are shown in Fig- 
ure [13] We find that the downsizing signal for M^ lv and 

M* piv is clearly detected in the COSMOS data using the 
methods of BIO. This provides an independent test on 
our detected evolution of the pivot quantities and sug- 
gests that the detected downsizing signal is robust to the 
methodology that is employed. 

This test also raises the interesting question of how the 
method used in this paper (an HOD-based model that 
includes fits to clustering and g-g lensing) compares to 
the method of BIO (abundance matching using only the 
SMF) . As can be seen in Figure [131 we find that the two 
methods yield very similar results for the pivot quanti- 
ties. We do however find subtle differences in the actual 
SHMRs between the two methods. Tracking down the 
cause of the exact differences, although a very interesting 
question in itself, is beyond the scope of this paper and 
we defer this study to follow-up work. For the purposes 
of this paper, we will simply emphasize that the analysis 
of BIO applied to the COSMOS results fully agrees with 
our claims concerning the evolution of the pivot quanti- 
ties. 

5.6. The role of the stellar mass function 

In our analysis, the errors on the SMF are small com- 
pared to the clustering and the lensing. It is always the 
case that a measurement of a one-point statistic from a 
given set of data is more precise than a measurement of a 
two-point (or higher) statistic. Thus, the SMF plays an 
important role in constraining our parameter set. There- 
fore, we investigate the SMF in further detail in this sec- 
tion, and in particular, we show a more in-depth com- 
parison with SDSS mass functions. 

Figure [14] shows the COSMOS mass functions com- 
pared to various SDSS mass functio ns that have been 
commonly employed in the literature (jPanter et al"]l2007t 
IBaldrv et alJ[200a ILi fe White! 120091) . The main differ- 
ence that we may expect between the COSMOS mass 
functions and the SDSS ones (besides sample variance 
and systematic error) is that the high end of the mass 
function may be inflated due to a larger value of <Ji og M. 
in COSMOS. To gauge how much the COSMOS mass 
functions arc affected by Eddington bias compared to 
SDSS, we use our model to predict the COSMOS mass 
functions, de-convolved to the expected scatter for SDSS 
(ciogM, ~ 0.17 dex). The results are shown in the right 
hand panel of Figure [14] We find that the difference 
in scatter is not significant enough to explain the differ- 
ences between the COSMOS and SDSS mass functions. 
It is more likely that the differences are due, for example, 
to varying assumptions regarding stellar population and 
dust models. 

The difference between COSMOS and ILi fe White! 
(|2009l ) corresponds roughly "by eye" to a "left/right" 
shift along the X-axis (log lp (M' i ) - log 10 (M™ smos ) - 
0.2). This difference is within the estimated systematic 
uncertainties (0.25 dex according to B10). However, this 
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type of systematic shift will be reflected directly in the 
SHMR (Figures [10] and El) by a "left /right" shift along 
the X-axis and will also affect the normalization of M* IV 
and (Mh/M*) piv . This ~0.2 dex shift would bring the 
normalization of BIO into closer agreement with our re- 
sults. Tracking down the exact source of this systematic 
shift is beyond the scope of this paper but it could be 
associated with differences in the assumed dust mode0 
for example. 

We conclude that in order to use the low-z SDSS data 
as a z ~ anchor to study the redshift evolution of the 
SHMR, a homogeneous analysis of both the SDSS and 
the COSMOS data is critical. This will be the focus of a 
future paper. 

5.7. The total galaxy stellar content as a function of 
halo mass 

Figure [T5l shows the conditional stellar mass functions, 
for various halo masses and redshifts from our best fit 
model. We can use these functions to calculate the total 
amount of stellar material locked up in galaxies as a func- 
tion of halo mass, noted hereafter M* ot (see Equation 16 
in Paper I). Investigating M' ot is of interest because it 
reveals the efficiency with which dark matter halos accu- 
mulate stellar mass from the combined effects of in-situ 
star formation and accretion via merging. 

One might worry that calculating the total amount of 
stellar material locked up in satellite galaxies requires ex- 
trapolating our model beyond the lower and upper stellar 
mass bounds for which our model has been calibrated. As 
discussed in § 15.5.11 the extrapolat ion of our model is in 
good agreement wit h results from | Blanton et al.l (|2008[ ) 
at low M* and with IHoekstral (|2007l) at high M*. Thus, 
to first order, this extrapolation does not appear unrea- 
sonable. Let us therefore make the assumption that the 
extrapolation of our SHMR is not wildly incorrect. We 
will now investigate which mass range of satellite galaxies 
contribute most to M* ot . 

From Figure [131 it is clear that satellite galaxies are 
a subdominant component of the total stellar mass at 
AI h = 10 12 M Q . Our present concern is therefore 
only relevant for M/, > 10 12 M Q . Let us consider 
the total stellar mass associated with satellite galax- 
ies as a function of Mh in a fixed stellar mass bin: 
M* ot ' sat (i\4|Mf ,Mf ). As shown in Paper I, the ex- 
pression for M* ot ' sat (M h |Mf,Mf) is given by: 

/•Mf 

Ml ot > s * t {M h \Ml 1 ,M?)= / $ fl (M*|M h )Af*dM*. 

(20) 

We have tested how M^'^M^M* 1 , Mf ) varies with 
the integral limits, M* 1 and M* 2 . We find that at fixed 
halo mass, most of the stellar mass associated with satel- 
lite galaxies arises from a relatively narrow range in stel- 
lar mass. In particular, for halos with Mh > 10 12 M Q , 
the bulk of M„° ' sa is built from satellite galaxies in the 
range 10 10 M Q < M* < 10 n Af Q . Therefore, provided 
that the extrapolation of our model is not wildly incor- 

8 Our stellar masses use d the lCharlot fc F all (2000) dust model 
whereas ITT fc White! 1120091 1 use lBlantorTg" Rowcis (2007). 



rcct, the bulk of M*° ' sa arises from satellites that are 
within the tested limits of our model. 

Having underlined this caveat, we have calculated M* ot 
using the best fit parameters for each of the three redshift 
bins and the results are shown in Figure [TBI This Figure 
will be discussed in detail in the following section. 

6. DISCUSSION 

Using a self-consistent framework to simultaneously fit 
the g-g lensing, spatial clustering, and number densities 
of galaxies in COSMOS, we have obtained a robust char- 
acterization of the evolving relationship between stellar 
mass and halo mass over two orders of magnitude in M* . 
The nature of this relationship, shown in Figure 1161 is 
not only a byproduct of cosmic mass assembly but is 
also shaped by the physical processes that drive galaxy 
formation, ultimately providing valuable constraints on 
both. In this section we begin by discussing various pro- 
cesses that shape the form of M^/Mh versus Mh- We will 
then introduce a simple framework for interpreting evo- 
lution in this relation by considering the relative growth 
of stellar mass as compared to the growth of dark mat- 
ter halos. Finally, we will discuss the observed evolution 
of the pivot quantities and we will show how a constant 
pivot ratio may imply that the mechanism responsible 
for the shut-down of star formation in massive galaxies 
may have a physical dependence on M^/Mh- 

6.1. The total stellar mass content of dark matter halos 

Figure [JJ)] separates the stellar content of the aver- 
age dark matter halo into a contribution from the cen- 
tral galaxy and a contribution from the sum of satel- 
lite galaxies. Central galaxies show a M*/Mh ratio that 
rises steeply to a maximum at Mh ~ 10 12 M© before 
decreasing somewhat more gradually in halos of higher 
mass. The fact that halos above this mass scale (at 
the redshifts considered) have cooling times longer than 
their dynamical times has been invoked by modelers for 
some time to help explain why cooling and star forma- 
tion shut down at the highest masses, with some re- 
finement d ue to the presence of so - called "cold - mode 
accretion" (IBirnboim fc Dekell [2001 IKeres et al.l l200l 
IBirnboim et al.l 120071: ICattaneo et al1l2006[ ). Wewill re- 
turn to the evolution of this mass scale at a later point 
in the discussion. 

Central galaxies strongly dominate the total stellar 
mass content at Mh < 2 x 10 13 M ("the central domi- 
nated regime"), including at the peak mass, Mh ~ 10 12 
Mq, while the stellar mass in satellites dominates at 
M h >2x 10 13 M Q ("the satellite dominated regime"). 

The transition between the two regimes is driven by 
the steep decline in M^ cn /M h at M h > 10 12 M . This 
decline occurs as the contribution from satellites begins 
to rise. One might then naturally ask if central galaxies 
in group-scale halos experience stunted growth simply 
because stellar mass is accumulating within the halo in 
the form of satellite galaxies, instead of merging onto the 
central galaxy. Figure Q1)] reveals that this is not the case. 
Indeed, the solid line in this Figure demonstrates that 
the total stellar mass fraction of halos declines at Mh > 
10 12 Mq. Thus, even if all satellite galaxies were allowed 
to rapidly coalesce at the center of the potential well, 
the central galaxies of group-scale halos would still have 
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Fig. 14. — Comparison of the COSMOS stellar mass functions to the SDSS mass functions of Pantcr ct al. ( 2007'). [Baldrv et al.l {2008), 
and lLi fc White] (200^ ). In the left panel we show the COSMOS z\ mass function. The error bars represent the expected sample variance 
in COSMOS. In the right panel we show the mass functions for all three COSMOS mass functions, de-convolved to a common o"i og M* °f 
0.17 dex. The differences between the COSMOS mass functions and the SDSS mass functions cannot be explained by Eddington bias and 
are therefore more likely to be due, for example, to differences in assumptions regarding stellar population and/or dust models. 



lower M*/Mh ratios than those in halos of Mh ~ 10 12 
M Q . Thus, we conclude that dark matter halos globally 
decline in the efficiency by which they accumulate stellar 
mass at M h > 10 12 M Q . 

We note that in massive halos, the intra-clustcr 
light (ICL) (not accounted for in this analysis) is esti- 
mated to con tribute an additional 20-30% of the total 
stellar mass (iFeldme er et al.l 120041: iZibetti et al.1 120051: 
iGonzalez et al.1 120051 : iKrick et al.1 120061 ). Adding this to 
the satellite component in Figure [TBI does little to bridge 
the factor of 2-4 gap in M* ot /Mh between the satellite 
component in high-mass halos and centrals at the peak 
mass. 

6.2. The role of galaxy mergers in determining the 
shape ofMl ot /M h 

The majority of the total mass in an average dark mat- 
ter halo is built f rom halo-halo merger s with mass ratios 
above 1:10 fe.g.. iHopkins et aLlETOlOal) . At halo masses 
below 10 12 M , the steep rise in Ml ot /M h with M h im- 
plies that the typical stellar mass ratio of galaxy mergers 
will be less than the typical mass ratio of the dark mat- 
ter halos hosting these galaxies. In other words, major 
halo mergers are minor galaxy mergers in this regime. 
Thus, the accumulation of stellar mass through the ef- 
fects of merging will be limited compared to the growth 
in total mass of such halos. The steep rise of Ml ot /Mh 
must therefore reflect the greater importance of star for- 
mation at masses bel ow Mh ~ 10* 2 over assembly 
from galaxy mergers (|Bundv et a l J 120091) . Similar con- 
clusions have also been reached bv iConrov fc Wechslerl 
see their Figures 2 and 3 in particular). 

Simple arguments suggest, however, that once halos 
grow past the pivot mass and in the absence of signifi- 
cant star formation, M* ot /M^ should dip below the peak 
value since these halos can only grow by merging with 
halos with lower values of M' ot /M^. At slightly higher 



mass, the decline in M* ot /M^, now means that stellar 
mass ratios are enhanced with respect to halo mass ra- 
tios, and the trend must reverse again. This repeating 
pattern should cause a flattening of M' ot /Mh above the 
pivot mass. While this behavior is certainly apparent in 
Figure [TBJ our different redshift bins also reveal that at 
fixed mass among high-mass halos (Mh > 4 x 10 13 ), the 
total stellar mass content declines at later epochs. We 
speculate that this trend could arise from the smooth ac- 
cretion of dark matter, which brings no new stellar mass, 
and amounts to as much as 40% of the growth of dark 
matter halos (e.g., Fakhouri & Ma 2010). One way to 
test this hypothesis would be to populate a z = 0.88 
N-body simulation with our 23 HOD and evolve the 
subhalo and halo populations to z = 0, assuming no 
star-formation. This would reveal the amount of stars 
that are acquired th rough mergers in this redshift range 
(jZentner et al.l [20051 ) . We note that the destruction of 
satellites and a growing ICL component could also con- 
tribute to this trend given that the ICL at z = could 
make up 20-30% of the stellar content of massive halos, 
roughly the amount by which M' ot /Mh declines over our 
redshift range. 

6.3. The pivot quantities and the quenching of star 
formation in central galaxies 

The location at which halos reach their maximum ac- 
cumulated stellar mass efficiency is encoded by the pivot 
mass quantities, Mf iv , Mf v , and (M h /M*) piv . While we 
observe downsizing trends for both M plv and M? lv , the 
co-evolution of these two parameters leaves (Mh/M*) plv 
constant with redshift. This can be seen in the right 
panel of Figure 1161 the pivot ratio evolves very little 
over our redshift range, while the mass scale of the peak 
(MP 1V ) does evolve downward by nearly a factor of 2. 
Given the low satellite content of halos at these masses, 
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Fig. 15. — Left panel: The conditional stellar mass functions for central and satellite galaxies for three different halos masses. Results 
are for z = 0.66. Solid curves represent <fr s (Af*|iW/j) while dashed curves represent <l> c (Af* | M^). The curves are normalized such that the 
total area is the mean number of galaxies at that halo mass. Right panel: Conditional stellar mass functions for halos of = 10 13 Mq 
for our three rcdshift bins. 



merging is not likely to play a dominant role in driving 
growth in M* below the pivot peak, arguing instead that 
the regulation of star formation is key to understanding 
this behavior. The physical process that sets the pivot 
peak at Mh = 10 12 M©, and drives the subsequent de- 
cline in M£ en /Mh at higher halo masses, must be linked 
to the shut-down of star formation in central galaxies. 
In the reminder of this discussion, we will focus on inter- 
preting the downsizing behaviour of the pivot quantities 
in this context. 

6.4. A simple model for interpreting evolution in 
M*/M h versus M h 

A complete and comprehensive interpretation of our re- 
sults requires modelling and accounting for dark matter 
accretion histories, galaxy merger rates, and star forma- 
tion r^te^_J^FR]_as_afunctifln of redshift (for example, 
see iConrov fc W cchsler 200!|). Nonetheless, we will in- 
troduce some toy models based on simple arguments to 
provide a first interpretation of our results. Our goal here 
is to evaluate our results in the context of other observa- 
tions and theoretical work on galaxy formation models 
and to set the stage for a more detailed treatment in 
subsequent work. We begin with a general treatment of 
evolution in the SHMR and then will focus on applying 
this treatment to interpret the evolution we observe in 
the pivot quantities. 

The physical basis for evolution in M*/Mh versus Mh 
must be considered carefully because the stellar mass 
and associated halo mass of a galaxy can evolve inde- 
pendently, depending on the mass scale involved and on 
processes including merging, smooth (diffuse) dark mat- 
ter accretion, star formation, and even tidal stripping. 
The sum of these processes on the growth of Mh and 
M* shapes the behavior of M*/Mh versus Mh in dif- 



ferent ways, as illustrated by the schematic diagram in 
Figure [17] Here we let -qi represent the ratio M*/Mh 
at redshift Zhi g h and rjo represent the equivalent ratio at 
rcdshift zi ow with Zhigh > z\ ovl . We further consider a 
dark matter halo of mass Mh that has grown by a rela- 
tive factor of \m h from Zhigh to z\ ow . We can write that 
M h (zi ow ) = X Mh x M h (zM S h)- Characterizing growth in 
the stellar mass of the central galaxy (although similar 
arguments apply to M* ot ) by a factor of Am,, we can 
simply write that: 



Am. = -xX Mh . (21) 

m 

If r/o > 771 , we infer that Aa/„ > Am (1 and that the stel- 
lar mass has experienced a stronger relative amount of 
growth compared to that of the dark matter from Zhigh 
to ziow If on the contrary, 770 < 771 then the relative 
growth of the stellar mass is less than the dark matter. 
Note that this schematic view demonstrates that identi- 
cal curves for M*/Mh versus Mh at different redshifts do 
not necessarily indicate a lack of evolution, since Am,, is 
always greater than zero. 

We now apply this intuitive framework to the evolution 
observed in Figure 1161 Considering values of 771 and 770 
applied to the total stellar mass curves in Figure HH we 
see that below Mh = 10 12 , the fractional growth in stellar 
mass outweighs the growth in halo mass. This reflects 
the greater importance of star formation at low masses 
over assembly from galaxy mergers, the same conclusion 
reached above by simply considering the shape of the 
SHMR. This evolutionary trend reverses above Mh = 
10 12 M©, consistent with the notion that star formation 
is largely shut down in centrals above this mass. 
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Fig. 16. — Total stellar content locked up in galaxies as a function of halo mass compared to the cosmic baryon fraction measured by 
the Wilkinson Microwave Probe (WMAP5; f b = = 0.171 ± 0.009: [Duiiklev et aI.ll2u09T ). Left panel: Our prediction from the z\ bin 
(z ~ 0.37). Right panel: Our three redshift bins. z\ is shown by the solid dark blue line, 22 is shown by the blue line, and 23 is shown 
by the turquoise line. Dotted lines show the contribution to M* ot from the central galaxy and dashed lines show the contribution from 
satellite galaxies. Shaded regions represent the errors on M" * ot /M2006 • M* ot is dominated by the central galaxy at M h < 2 X 10 13 Mq and 
by satellites at M h > 2 X 10 13 Mq. 
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Fig. 17. — Schematic illustration of our results in terms of galaxy 
mass assembly versus halo growth. The turquoise line represents 
our 23 result and the dark blue line represents our 21 result. From 
the relative positions of these two curves we can infer that from 
2 ~ 0.88 to 2 ~ 0.37 and at M h < 10 12 M Q , the stellar mass of 
the central galaxy has experienced a stronger growth in proportion 
to the growth of the dark matter. On the contrary, at M h > 10 12 
Mq , the stellar mass of the central galaxy has experienced a more 
mild growth in proportion to the dark matter. 

6.5. Understanding the evolution of the pivot quantities 

We now apply these simple arguments to the evolution 
in the pivot quantities. We focus only on central galax- 
ies, neglecting the minor contribution from satellites near 
the pivot mass. The aim here is to explore several sim- 
ple models for the quenching of star formation and to 
investigate which models might reproduce the observed 
evolution of the pivot quantities, namely, a pivot halo 
and stellar mass that decrease at later epochs (downsiz- 
ing) but leave the pivot ratio constant. 



Shown schematically in Figure [TBI we consider how 
our high redshift (z = 0.88) SHMR relation would 
evolve towards lower redshifts under several prescrip- 
tions for stellar and halo growth. We begin with no 
assumptions about the SFR but adopt a halo growth 
rate (\M h ) that is roughly constant over the mass range 
spanned by the peak. This assumption is well justified 
by dark matter m ass accretion rates derived from N- 
body simulations ( Wechsl cr et al.l I2002t IMcBride et al 



I2009t iFakhouri et al.ll2010D . For example , IFakhouri eta! 



(|2010D find that M^jMh is only weakly dependent on 

halo mass with M h /M h cx Mfc 1 . 

We now consider several different quenching models 
and investigate their impact on the redshift evolution of 
the pivot quantities. 

• No quenching model: to begin with, we consider a 
model with no quenching of star formation and in 
which the stellar growth rate (Am, ) is constant over 
the halo mass range spanned by the peak (Row A 
in Figure HH]). We consider the redshift evolution 
of the SHMR for three values of Am, defined with 
respect to Am,, • In all three cases, this model leads 
to an increase in M£ IV with time (contrary to what 
we observe). We can therefore conclude that Am, 
must vary with Mh, not a surprise given the expec- 
tation that the SFR shuts down above M/, = 10 12 
M . 

• Fixed halo mass for quenching: we next consider 
a model in Row B in which star formation is 
quenched as galaxies cross above a fixed halo mass, 
M q . However, this instantaneous quenching model 
fails to reproduce (not surprisingly) the downward 
evolution of Mj^ lv and also yields evolution in the 
pivot ratio, which is not detected. 
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Row A 

Simple evolution 
model: a constant 
stellar and halo 
growth rate across 
the peak. In this 
case, the pivot 
mass always 
increases at later 




times. lo 9( halo mass ) log( halo mass ) log( halo mass ) 
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Fixed halo mass „ 

(Mq) above which J 

star formation is -5 

quenched. For 2 

Mh>Mq, XM*=0. a 

This scenario fails — 
to produce downsizing. 



Row C 

Redshift dependent 
Mq above which star ^ 
formation is j 
quenched. The 
middle panel 2 
produces downsizing ~5i 
in the pivot halo — 
mass, but does not 
necessarily keep the 
pivot ratio constant. 
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Row D 

Schematic version 
of our results. We find 
downsizing in the pivot 
stellar and halo mass and 
a constant pivot ratio. 
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Row E 

Model in which star 
formation is limited 
by a critical pivot 
ratio, Hcrit- In this 
model, the middle 
panel reproduces 
our results. 
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Fig. 18. — A schematic (and simplistic picture) of how Mt/M^ (for central galaxies) varies with redshift for different quenching models 
and for various prescriptions for stellar growth and halo growth (parametrized here by Am* and ^M h respectively). Our 23 (z=0.88) relation 
is represented by the solid red line. In this picture, we are not interested in understanding why the high redshift relation has the particular 
form that is observed, but simply in predicting roughly what the evolution of this relation should look like given various quenching models. 
The orange dashed line and the blue dotted line show how we expect the SHMR to evolve with time. In Row A we consider a model 
without any quenching and where Am* an d ^M h are constant with redshift and halo mass. In this model, the pivot halo mass will increase 
at later epochs which is not what we observe. In Row B, we consider a model in which star formation is quenched above a fixed halo 
mass, M q . This model fails to produce the observed downsizing behaviour for M? . In Row C, we consider a model where M q decreases 
at lower redshifts. In this case when Am* > ^M h below M q (Row C, middle panel), M^ lv follows a downsizing trend. However, in this 
scenario, a pivot ratio that is constant with redshift requires a fine tuning between Am* and the rate at which M q declines. In Row E, 
instead of assuming that star formation is quenched at a fixed halo mass, we now assume that star formation is quenched at a fixed critical 
Mt/Mh ratio (rj cr i t ). Given this assumption and if Am* > ^M h below f7 cr it (a reasonable assumption given estimates for halo growth and 
star formation at these scales), we can qualitatively reproduce our main results (compare the middle panel of Row D to the middel panel 
of Row E). 
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• Rcdshift dependent halo mass for quenching: Row 
C shows a model in which M q shifts downward with 
time. This model leads to downsizing in the pivot 
halo mass if Am, > Am,, below M q (Row C, middle 
panel). However, in order to keep the pivot ratio 
fixed in this scenario, the growth rate, Am, , must 
be tuned with respect to the rate at which M q de- 
clines. This would require a fortuitous coincidence, 
but obviously cannot be dismissed as an explana- 
tion. 

• Critical M^/Mh ratio for quenching: we finally 
consider an alternative scenario in which star for- 
mation is limited by a critical mass ratio, 7y cr it = 
M*/M h ps 0.04 (Row E). In this model, we can 
qualitatively reproduce our main results, including 
the downsizing trends in the pivot stellar and halo 
mass, and, by construction, a constant pivot ra- 
tio set by rfcrit (compare the middle panel of Row 
D to the middle panel of Row E). However, in or- 
der to produce downsizing behaviour in this model, 
Am, must be larger than Am ;i below r? cr it. It is im- 
portant to note that if Am, < Am,, below r} CT it, this 
model would fail to produce downsizing. 

The model explored in Row E seems a promising and 
simple mechanism that can explain the observed evolu- 
tion of the pivot quantities. We now test if observations 
are consistent with the requirement that Am, > ^M h 
below M h ~ 10 12 M Q . A halo of mass M h ~ 2 x 
IQ^Mq gr ows by a factor of Ay h ~ 1.4 from z = 0.88 
to z = 0.37 (|Fakhouri et al.ll2010f ). Star formation rates 
as a funct ion of Af» and redsh ift have recently be e n mea - 
sure d bv iNoeske et all j2007h . iCowie fe Barged IpOOl . 
and iGilbank et al.l (|2010T l . These measurements indicate 
that galaxies of mass M* ~ 10 10 M Q grow by roughly a 
factor of Am, = 2 - 3 from z = 0.88 to z = 0.37. There- 
fore, at z < 1 and for Mh < M^ lv , current estimates for 
halo growth coupled with estimates for stellar growth are 
indeed consistent with our observation that stellar mass 
has experienced a stronger relative amount of growth 
compared to that of the dark matter (Am,, < Am, be- 
low 77crit)- 

We further note that the relatively weak dependence 
of Am, on Mh required to maintain a constant slope in 

M*/Mf, versus M h with redshift (again for M h < M£ iv ) 
is implied by the weak SSFR-A-/* relation observed for 
galaxies with M* < M* lv (a typical power-law fit gives 
SSFR- Af° 4 ). This could explain why (3 (the low mass 
slope of the SHMR) is observed to remain remarkably 
constant at /3 = 0.46 at z < 1. 

Given the simple arguments outlined above, we argue 
that the fact that the pivot ratio remains constant may 
suggest that star formation is fundamentally limited by 
a critical mass ratio, r/ cl i t = M*/M^ ps 0.04. A very el- 
egant and compelling consequence of this model is that 
the observed downsizing trends in M* lv and M^ lv can 
be automatically explained given the observation that 
Am, > A Mi, below Mh ~ 10 12 M©. Previous work by 
iBundv et al.l (|2006[) on the evolution of the galaxy stel- 
lar mass function has found a similar downsizing trend 
for the "transition mass", M tr , which is defined in terms 
of M* by the declining fraction of star-forming galax- 



ies at the h i ghest masses. The transition mass from 
IBundv et all (|2006[) evolves from M tr - 9 x 10 10 M 
from z ~ 0.9 to M tr - 5 x 10 10 M from z ~ 0.6. 
This transition mass is similar to M|" v , reinforcing the 
notion that the pivot mass marks the end of rapid star 
formation among h alos. The obvious difference between 
IBundv et al.l (|2006|) and this paper is that our work adds 
a key missing ingredient, which is the evolution of the 
pivot halo mass and the pivot ratio. If the quenching 
of star formation depends on a critical M*/M/j ratio 
then the fact that low mass galaxies grow more rapidly 
than dark matter below the pivot scale provides a simple 
explanation for the observed downsizing i n the sites of 
star fo rmation observed by studies such as IBundv et al.l 
(l200l . 

6.6. Physical mechanisms that might depend on M^/Mh 

In the previous section, we demonstrated that a con- 
stant pivot ratio provides important clues concerning the 
physical mechanisms that quench star formation. We 
now discuss possible mechanisms that might tie quench- 
ing to M*/M h . 

The notion of a fixed maximum stellar-to-dark matter 
ratio, ry cr it, has been relatively unexplored in the liter- 
ature. Theoretical arguments tend t o favor a relatively 
fixed (if broad) critical halo mass (see iBirnboim fc Dekell 
120031) at z < 1, with significant modifications from 
cold-mode accretion occurring mostly at higher redshifts. 
But, quenching at fixed halo mass alone is not suffi- 
cient to reproduce the local stellar mass function and 
red-sequence fraction, and also fails from simple argu- 
ments to produce downsizing in the pivot masses, as 
shown by Row B of Figure [T5] As a result, most semi- 
analytic models include a quenching channel initiated by 
mergers or disk instabilities among galaxies in halos be- 
low the critical halo ma s s threshold (e.g.. iBower et al.l 
12001 iCroton et aLll200l iCattaneo et all 120061) . Above 
the critical halo mass, once a quasi-static halo of hot 
gas has formed, low-luminosity feedback (i.e., radio- mode 
AGN feedback) is often invoked to prevent cooling and 
star-formation in massive halos at late times. Using sim- 
ple arguments, we have shown that mergers do not pro- 
vide a likely explanation for determining the scale of the 
pivot mass. We therefore focus on AGN feedback and 
disk instabilities as possible mechanisms. 

It is common practice to consider low-luminosity AGN 
feedback only in halos above a fixed halo mass. However, 
in practice, efficient AGN feedback is more complex and 
requires at least two ingredients. First, a quasi-static 
halo of hot gas must exist to which the AGN jets can 
couple. However, AGN feedback also requires a suffi- 
ciently large black-hole to produce jets powerful enough 
to initiate this coupling. Given that gas cooling rates 
scale roughly with halo mass and that black hole mass 
scale s roughly with galaxy mass, via the — a r ela- 
tion (jGebhardt et alj|2000t IPerrarese fc Merrittll2000[) . it 
is reasonable to speculate that AGN feedback efficiency 
would depend on M*/M/j. 

Considering the galaxy population broadly, this sce- 
nario requires a sufficiently large bulge component for 
AGN quenching to be effective. While bulges may be 
built stochastically in galaxy mergers, a further link may 
tie secular bulge formation via disk instabilities to the 
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value of M*/Mh, thereby cementing the relationship be- 
tween quenching and Tfcrit- It has been shown that disks 
become unstable to bar modes i f the disk mass domi- 
nates the gravi tational potential (jEfstathiou et alJl!982t 
IMo et aI]|1998D . Semi analytic models typically consider 
that disk instabilities occurs if: 

V max /(G N M disk /r disk ) - 5 < 1, (22) 

where Mdisk represents disk mass, rdisk is the disk ra- 
dius, and Vmax is the maximum of the rotation curve. 
Depending on the implementation of this criterion, Vmax 
may be equated either to the halo virial velocity, or the 
disk velocity at its half mass radius (see discussion in 
iParrv et al.Tl2009T ). An instability will cause either a par- 
tial or a total collapse of the disk, leading to a burst of 
star formation at the center, the formation of a spheroid, 
and also possibly fueling the central black hole. The disk 
instability criterion in Equation [22] shows a dependence 
on (Af/j/Mdisk) (relating V max to Mh) which could be 
reflected in M^/M* as the instability converts cold gas 
into stars. Disk instabilities might therefore play a role 
in setting the pivot masses and enforcing r/ cr i t . The co- 
incident fueling of the central black hole during the in- 
stability may help initiate AGN quenching and regulate 
the global decline in M^ en /Mh beyond the pivot mass. 

Assuming the disk instability framework, we can test 
whether the predicted sizes of stable galactic disks given 
our critical pivot ratio, 7y cr it = M^/Mh « 0.04, are con- 
sistent with observations. Approximating the halo virial 
velocity for Vmax, we rewrite Equation [22] to derive the 
maximum size of stable disks: 

r dis k < Rh x (M*/M h ). (23) 

Applying this criteria to the pivot masses yields the 
condition that rdisk < 8 kpc. Interestingly, this con- 
dition is very well satisfied by the observed size distri- 
butions of disk galaxies which tend to fall off rapidly 
just below rdisk = 8 kpc at both z ~ (see Figure 
11 inlShen et all 120031) and at z > (see Figure 10 in 
iSargent et al.ll2007l ). We conclude that, via the growth 
of bulges and initialization of AGN feedback, disk insta- 
bilities provide a promising link between our observed 
critical pivot ratio and quenching, although further in- 
vestigation is clearly needed. 

Finally, we note that ?y cr it might also be related to the 
competition between the cooling and accretion of cold gas 
in central star-bursts and the resulting feedback from ei- 
ther star formation itself or a co-evolving quasar-mode 
AGN. This Eddington-likc li mit has been explored i n the 
context of stellar systems by iHopkins et al.l (]2010bf ) who 
derive a maximum stellar surface density from simple ar- 
guments. In the context of dark matter halos explored 
here, similar arguments might naturally yield a fixed 
value for r) CI - lt similar to that obtained by our analysis. 

7. CONCLUSIONS 

The aim of this paper is to study the form and evo- 
lution of the stellar-to-halo mass relation (SHMR) from 
z = 0.2 to z = 1.0. To achieve this goal, we have per- 
formed a joint analysis of galaxy-galaxy lensing, spatial 
clustering, and number densities of galaxies in COSMOS. 
As a result, we have obtained a robust characterization of 
the evolving relationship between stellar mass and halo 



mass over two to three orders of magnitude in M*. The 
nature of this relationship is not only a byproduct of 
cosmic mass assembly but is also shaped by the physical 
processes that drive galaxy formation, ultimately pro- 
viding valuable constraints on both. A complete and 
comprehensive interpretation of our results requires mod- 
elling and accounting for dark matter accretion histories, 
galaxy merger rates, and star formation rates as a func- 
tion of redshift. This will be the focus of future work. 
Nonetheless, we show how simple evolutionary models 
can already provide a first interpretation of our results, 
setting the stage for a more detailed treatment in future 
work. Using simple arguments, we evaluate our results 
in the context of other observations and theoretical work 
on galaxy formation models. 

We have defined the pivot quantities (Mf IV , M piv , and 
(M? l /M») plv ) as the location at which halos reach their 
maximum accumulated stellar mass efficiency. The evo- 
lution of the pivot quantities contain key clues about the 
physical processes that are responsible for the quench- 
ing of star formation in halos above Mh > 10 12 M Q . 

While we observe downsizing trends for both M plv and 
M^ lv , the co-evolution of these two parameters leaves 
(M/ l /M*) plv roughly constant with redshift. We argue 
that this result raises the intriguing possibility that the 
quenching of star formation may have a physical depen- 
dence on Mh/M* and not simply on M% as is commonly 
assumed. If the quenching of star formation indeed de- 
pends on a critical M*/Mh ratio then the fact that low 
mass galaxies grow more rapidly than dark matter be- 
low the pivot scale provides a simple explanation for ob- 
servations of downsizing in the sites of star formation. 
Additional and more precise measurements of the pivot 
quantities would be highly interesting in order to con- 
firm whether or not the pivot ratio remains constant, or 
if instead it evolves mildly with redshift. Interestingly , 
there are hints from the results of iBehroozi et all (j2010f ) 
that the pivot mass might in fact remain constant back 
to z = 4. 

Wc highlight four avenues of exploration that would 
be interesting to pursue in the future and that would 
improve this analysis. Firstly, we note that the compari- 
son of our results with SP SS abundance matching results 
(e.g. IBehroozi et aTll2010D is limited by systematic differ- 
ences between stellar mass estimates. In order to use the 
SDSS data-set as a low-z anchor point, a homogeneous 
analysis of SDSS and COSMOS would be necessary. Sec- 
ondly, our work in COSMOS is limited at the high mass 
end by sample variance. Larger data-sets at higher red- 
shifts than SDSS such as the Baryon Oscillation Spectro- 
scopic Survey (BOSS0 and the Canada-Francc-Hawaii 
Telescope Legacy Survey (CFHTLS0 should provide in- 
teresting constraints on the evolution of the high mass 
end of the SHMR. Thirdly, pushing the SHMR down 
to 10 8 M© < M* < 10 9 M using techniques such as 
described in this paper is clearly an exciting avenue to 
explore and will be facilitated by upcoming very deep 
optical and IR surveys such at UltraVista and the Hyper 
Suprime Cam (HSC) survey on the Subaru telescope. Fi- 
nally, we note that the stellar mass function is a powerful 

9 http://cosmology.lbl.gov/BOSS/ 

10 http:/ /www. cfht.hawaii.edu/Scicnce/CFHLS/ 
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tool for placing constraints on the SHMR. However, the 
derivation of SMFs is clearly currently limited by sys- 
tematic uncertainties in stellar mass estimates. Working 
towards an improved understanding of stellar mass es- 
timates and towards reducing systematic errors in the 
SMF will be the single most important avenue for im- 
proving the type of analysis presented in this paper. 

Finally, while our analysis demonstrates that the com- 
bination of multiple and complementary dark matter 
probes is a powerful tool with which to elucidate the 
galaxy-dark matter connection, we emphasize that such 
probe combinations also hold great potential to constrain 
fundamental physics, including the cosmological model 
and the nature of gravity. Exploring the sensitivity of 
the combination of g-g lensing, clustering, and the SMF 
to cosmological parameters will be the focus of a follow- 
up paper. 
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